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CHAPTER I 


INTRODUCTION 

The increasing cost of fuel is exerting great pressure 
on aircraft manufacturers to develop more energy efficient 
aircraft. Many of the refinements and improvements neces- 
sary in achieving more efficient aircraft will be by appli- 
cations of Active Controls Technology (ACT) . In this con- 
cept aircraft efficiency improvements are made by the in- 
creased dependence on augmentation control systems. 

The largest benefit of active controls will occur when 
the design process incorporates ACT concepts into airplane 
configurations at the preliminary design stage. However, 
this change in design philosophy is not likely to occur un- 
til much more reliability is achieved in predicting aero- 
dynamic and structural characteristics accurately. The air- 
craft industry is relying on government and military funded 
research to advance the technology needed to enable reliable 
predictions of aerodynamic, structural, and control system 
behavior. 

Of all ACT concepts , flutter suppression has the high- 
est technical risk. The consequences of a design error 
from a safety standpoint are catastrophic. The National 
Aeronautics and Space Administration (NASA) is actively 
pursuing applications of active controls in a technology 
development program that will provide the data which will 



lead to a better understanding of unsteady aerodynamics and 
structural dynamics required to alleviate some of the risk 
in flutter suppression designs. Plight test evaluations of 
active flutter suppression controls will be performed in a 
program called DAST (Drones for Aeronautic and Structural 
Testing) which utilizes a highly modified BQM-34E/F Firebee 
II supersonic target drone as a testbed. 

Flutter Suppression Systems (FSS) have been demonstrat- 
ed on exisiting aircraft such as the B-52 CCV (reference 1) 
but these have had known and well defined flutter modes. 

What makes the DAST program unique is the fact that the FSS 
is designed to suppress a flutter mode whose characteristics 
have only been predicted analytically. 

The objectives of the flight test activity for the 
flutter suppression experiment, known as DAST I, are: 

• Validate active controls technology for 
flutter suppression, 

• Enhance and verify transonic flutter predict- 
ion techniques. 

A secondary objective is to provide a data base for tran- 
sonic aerodynamic loads prediction techniques for elastic 
structures. 

The validation of ACT for flutter suppression will be 
accomplished by demonstrating in flight a 20 per cent in- 
crease in the flutter speed of a transonic flutter critical 
vling and correlate the results with theoretical data. 
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Only in the relatively recent past have attempts been 
made to analytically predict flutter characteristics in the 
transonic speed range. The flight test activity of the DAST 
program will provide the data needed for evaluation of the 
developing transonic flutter prediction techniques. 

This report will describe the development of the DAST 
I aircraft and some of the unique aspects of testing this vehicle. 

Chapter 2 will provide a description of the DAST I 
remotely piloted research vehicle (RPRV) and the support 
requirements for flight operations. In Chapter 3 the pre- 
dicted flutter and rigid body flight characteristics are 
presented. The synthesis of the flutter suppression con- 
trol system, the primary flight control system, and the 
back-up autopilot are described in Chapter 4. Chapter 5 
describes the hardware and software implementation of the 
DAST I vehicle systems. 

The preparation for the flight testing phase of the 
DAST program is presented in Chapter 6, This includes a 
discussion of testing philosophy and methods. A simulation 
of the flexible and rigid body modes of the DAST I aircraft 
is used to evaluate the near-real time testing procedures 
that will be used for the FSS flight experiment. 

The conclusions of this report are presented in 
Chapter 7 along with recommendations for further study. 
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CHAPTER 2 


OPERATIONAL APPROACH 

A description of the DAST I vehicle and its support 
requirements are presented in this chapter. The vehicle's 
estimated operating envelope is defined and a description 
of the remotely augmented control mode is given. The basic 
elements involved in flight testing this aircraft are pre- 
sented. 

2.1 Test Vehicle Description - The flight testing 
phase of the DAST I program will involve a highly modified 
BQM-34E/F supersonic target drone shown in figure 2.1. This 
vehicle has a supercritical wing (SCW) , designed to flutter 
within the operating envelope of the aircraft, in place of 
the standard drone wing. The wing, called ARW-1 (Aeroelas- 
tic Research Wing - 1) , is a 3/8 scale of the F-8 SCW (ref- 
erence 2) which is representative of possible transport air- 
craft wings designed to cruise at near sonic velocity. 

Small control surfaces on the trailing edge of the wing near 
the tip are utilized in the active flutter suppression sys- 
tem. 

No external modification to the DAST fuselage was re- 
quired, however, extensive internal changes for the new 
avionics and other systems have been made. The internal 
system changes are listed below: 

• Addition of FSS electronic and hydraulic 
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FIGURE 2.1 BQM-34E/F DRONE THREE VIEW WITH SUPERCRITICAL WING (DRST I) 


systems . 

• Replacement of the wing tip antennas with 
wrap-around nose cone antenna. 

• Replacement of the standard drone discrete 
autopilot with a simplified analog proportion- 
al system. 

• The addition of avionics for primary digital 
proportional remotely augmented control system. 

• Replacement of standard vehicle sensors with 
new instruments and PCM systems. 

Further details on these hardware changes are presented in 
Chapter 5. 

The DAST I estimated flight envelope is seen in figure 
2.2. This performance envelope is based on the engine cap- 
abilities defined in reference 3 and drag characteristics 
estimated from wind tunnel tests of a 1/6 scale DAST I mod- 
el. Drag information at Mach numbers greater than 0.98 
are based on data from the F-8 SCW aircraft (reference 2) . 

The 80 per cent throttle operating minimum seen in figure 2.2 
is a placard based on safe continuous engine operating con- 
ditions. An unsteady pitch-up occuring at high angles of 
attack will occur above the placard also indicated in figure 
2.2. Although the DAST I has the capability to fly super- 
sonically, no plans exist to fly the aircraft at speeds 
greater than M = 0.98. 
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The vehicle can carry 138.32 kilograms of fuel which 
gives it a flight time duration of approximately one hour 
at high altitudes. The fuel consumption rates increase 
dramatically at altitudes below 7,600 meters, however, and 
flight time duration at 3048 meters and M = 0.98 is only ap- 
proximately 13 minutes. 

2.2 Description of RPRV Operations - The support re- 
quirements to operate the DAST 1 as a RPRV at NASA DFRC are 
listed below: 

• RPRV launch aircraft 

• RPRV ground facility 

• Secondary command station 

• RPRV recovery vehicle 

Each of these elements are described in Sections 2.2.1 through 

2.2.4. 

The advantages of the RPRV testing technique are de- 
monstrated in high risk experiments such as the DAST I pro- 
gram. The development costs of full scale manned aircraft 
are eliminated by building simplified subscale RPRV's 
which if destroyed, neither risk the pilot nor result in a 
large financial loss. Public safety is not jeopardized be- 
cause of the remote testing site at NASA DFRC in the Mojave 

Desert in California. Flexibility in critical flight sys- 
tems is achieved by utilizing the facilities already develop- 
ed in the remote augmented vehicle (RAV) concept. 
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2.2.1 RAV Concept - The primary flight control system 
for the DAST I will be by remote augmentation. In this con- 
cept the control laws which augment the flying characteris- 
tics of the aircraft are programmed in FORTRAN in a general 
purpose computer on the ground. Closed loop operation is 
achieved by a telemetry uplink/downlink between the vehicle 
and the ground facility. This operation is illustrated in 
figure 2.3. The pilot has a full instrument display in the 
ground cockpit and plot boards indicating the vehicle's 
flight pattern and altitude are available to indicate de- 
viations in the flight plan. The details of the mechaniza- 
tion of the primary control system are given in sections 

5 . 4 and 5.5. 

The RAV concept has been successfully demonstrated on 
a 3/8 F-15 RPRV (reference 4) and as part of the F-8 Digit- 
al Fly-By-Wire aircraft presently being flown at NASA DFRC. 

2.2.2 F-104 Secondary Command Station - Air Force re- 
quirements at the Edwards, California test range specify 
that all powered RPRV's need a secondary and independent 
command station in the event of a loss in the telemetry 
link between the ground station and the test vehicle. For 
the DAST program this function will be performed by a 2- 
seat F-104 chase aircraft. In the rear seat of the F-104 

a small control panel is installed from which the chase air- 
craft can command the DAST in much the same manner as the 
ground facility. However, only an uplink signal to the 
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FIGURE 2.3 REMOTELY AUGMENTED VEHICLE CONCEPT 








DAST is active. Stability augmentation is supplied by an 
on-board autopilot in the RPRV. 

Figure 2.4 shows the range capability of the F-104 
transmitter. Oil is ejected into the jet exhaust of the 
DAST to enhance visable tracking of the vehicle in the oc- 
curance that the chase plane loses sight of the DAST vehi- 
cle . The secondary command station will fly the RPRV to 
a recovery area since the loss of the telemetry command 
link terminates reception of the test data also. 

2.2.3 B-52 Launch Aircraft - The basic drone was de- 

signed to be air launched from an aircraft or ground launch- 
ed using jato rockets. The ARW-1 was not designed to with- 
stand the stresses of ground launch so the DAST 1 will be 
air launched from a Boeing B-52 aircraft. This B-52 air 
launch capability was demonstrated with a basic Firebee ve- 
hicle. The basic Firebee vehicle is shown in figure 2.5 
on the B-52 pylon just prior to launch. 

A launch panel operator in the B-52 monitors the RPRV’s 
systems and starts the drone's engine before the B-52 climbs 
above 3000 meters. When the RPRV has satisfied system 
checks performed by the ground facility, F-104 chase air- 
craft, and launch panel operator, the B-52 flies to the de- 
signated drop zone for launch. Before the vehicle is re- 
leased from the pylon, the RPRV's fuel tanks are "topped 
off" to insure that the vehicle is released with full tanks. 
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The RPRV is launched with its tail surfaces locked in 
a nose down command for the first three seconds after re- 
lease. This precaution is to guarantee adequate separation 
from the B-52. 

2.2.4 Mid-Air Recovery - The DAST I flights will ter- 
minate in a mid-air recovery (MARS) of the vehicle by heli- 
copter. At the end of a flight or if an emergency occurs 
the vehicle's recovery mode is activated. This starts the 
sequence listed below: 

1. Vehicle is transferred to autopilot control. 

2. Engine shutdown. 

3. The vehicle begins a moderate climb. 

4. Drogue chute deployment and systems shutdown. 

5. Main chute deployment. 

6. Helicopter capture of the drogue chute tow 

line. 

This procedure is illustrated in figure 2.6. Events 
1 and 3 serve the purpose of starting internal timers and 
slowing down the aircraft. If speed requirements have been 
met (M <• 94) at the expiration of the drogue chute timer, 
the drogue chute is deployed and the vehicle's systems are 
turned off. Main chute deployment occurs when a second timer 
expires and the vehicle descends below 4570 meters. 

If a safety timer expires and events 4 or 5 have not 
taken place, the command to deploy the chutes is given re- 
gardless of flight conditions. Air Force experience with 
MARS technique has demonstrated an 80 percent success rate. 
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FIGURE 2.6 MID-AIR RECOVERY SEQUENCE 
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CHAPTER 3 


PREDICTED UNAUGMENTED CHARACTERISTICS 

This chapter presents the information defining the pre- 
dicted flutter and rigid body characteristics of the DAST I 
vehicle. These estimates have increased signifiance in the 
DAST flight test program since the active controls and 
flight controls were designed for these predicted character- 
istics. The planning for the flight tests will be done und- 
er the assumption that the estimated characteristics are 
accurate enough to preclude the requirement for further de- 
finition . 

3.1 ARW-1 Flutter Characteristics - Flutter is an 
aeroelastic phenomena characterized by a dynamic instability 
occuring in flight when the elasticity of the structure 
plays an essential part in the instability. The DAST I 
wing has been designed so that flutter will occur within 
its operating envelope. The problem of predicting the 
flutter characteristics and designing the FSS involves an 
accurate determination of the ARW-1 structural and unsteady 
aerodynamic characteristics. 

The development of the analytical model of the un- 
steady airloads and equations of motion are described in 
section 4.1.1. The representation of the DAST I structure 
used in the analysis is described below. 

3.1.1 Structural Dyanmics - The DAST I wing structure 
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consists of two machined steel spars, formed aluminum ribs 
spaced 30.84 centimeters apart, and fiberglass skins. A 
more detailed description of the wing is given in section 
5. 1. 

The vibration model for this wing that was used in the 
analysis was a NASA supplied NASTRAN model of the wing and 
wing center section structure. The ARW-1 idealization in- 
cluded: 

1. The removeable leading and trailing edge 
fiberglass structure, constructed of elements 
providing stiffness only for translational 
degrees of freedom. 

2. The ARW-1 skins were modeled as shear ele- 
ments with axial rods added to supply the 
membrane stiffness. 

3. The wing center section was modeled as hori- 
zontal beam and plate elements. 

4. The fuselage and empennage were modeled as 
beam elements with body mass represented as 
lumped masses. 

An illustration of this DAST structural model is 
shown in figure 3.1. The NASTRAN data for ARW-1 is found 
in reference 5. However, several changes were required to 
this NASTRAN model to more accurately describe the actual 
vehicle. The fuselage attachment of the wing was improved 
and places where the wing was very thin, foam core was 
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added to provide stiffness. A summary of improvements made 
in the NASTRAN model follows: 

1. Modification of the multi-point constraint 
equations to make the side of the fuselage 
rigid in pitch. 

2. Extension of spar webs past the solid alumi- 
num rib near the tip to model the vertical 
stiffness provided by the fiberglass tip's 
foam core. 

3. Addition of aluminum triangular membrane ele- 
ments at the trailing edge of the ribs to mod- 
el the foam core of the fiberglass trailing 
edge. 

4. The elastic modulus and shear modulus of the 
wing skin fiberglass material was found to be 
considerably different than had been estimated. 
Tests showed that it was necessary to decrease 
the elastic modulus by 14 percent and increase 
the shear modulus by 26 percent. 

5. A 0.907 Kg mass ballast was added approximate- 
ly 8.4 centimeters aft of the wing tip rib to 
compensate for the revised fiberglass charac- 
teristics and insure flutter within the flight 
envelope. 

Vibration solutions for three rigid body modes and ten 
elastic modes were found for both symmetric and antisymmetric 
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conditions. The predicted vibration modes are presented in 
Table 3.1 and 3.2. These values are compared to the pre- 
liminary results of a ground vibration test (GVT) conducted 
at NASA LaRC. Reasonable agreement exists between the ARW-1 
modes but considerable discrepancies occurred in the fuse- 
lage modes. The fuselage used in the GVT was that of the 
basic drone with standard equipment. The measured fuselage 
modes were also in considerable disagreement with GVT tests 
done on numerous standard drones for the military. For this 
reason, no modifications to the fuselage structural model 
will be done until further testing determines the source of 
discrepancies in the LaRC GVT. 

3.1.2 Flutter Boundary - The flutter mechanism of 
ARW-1 is a coalesence of the first and second wing modes. 
Both symmetric and anti- symmetric flutter modes exist with 
ARW-1. Root contours showing the two wing modes involved 
in the flutter are shown in figures 3.2 and 3.3. The sym- 
metric flutter mode locus also shows a vertical bending 
fuselage mode at approximately 20 hertz. It was included 
in the flutter analysis due to its correspondence to the 
flutter frequency. Figure 3.2 shows that the fuselage mode 
is slightly de-stabilized but does not become divergent be- 
fore the wing flutter mode. 

The ARW-1 flutter characteristics are described as 
explosive in nature. This is evident in figure 3.2 when 
the de-stabilizing movement of the wing roots are compared 
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TABLE 3.1 MODE DESCRIPTION - SYMMETRIC 
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FIGURE 3-4 DRST ARW-1 PREDICTED FLUTTER BOUNDARY 
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to changes in Mach number . Small changes in Mach number re- 
sult in large changes in the real and imaginary components 
of the elastic modes. Both the symmetric and anti-symmetric 
flutter characteristics are explosive in nature with the 
symmetric flutter occuring first at altitudes above 4270 
meters and the anti-symmetric mode fluttering first below 
4270 meters. This can be observed in figure 3.4. 

It is interesting to note from figures 3.2 and 3.3 
that at lower altitudes the second wing mode is the one 
which goes unstable while at the higher altitudes the wing 
bending mode encounters the instability. No instabilities 
in the higher frequency modes were predicted from the anal- 
ysis . 

3.2 PAST I Rigid Body Characteristics - The modifica- 
tion of the basic Firebee II flight simulator at NASA DFRC 
to account for the ARW-1 aerodynamics is described in this 
section. The digital six degree of freedom simulation re- 
quired rigid body aerodynamic coefficients in a format de- 
fined by angle of attack and Mach number break points. This 
format allows for non-linear characteristics to be included. 
The data described in section 3.2.1 was used in the flight 
planning simulator and the flight controls system analysis 
of sections 4.2 and 4.3. 

3.2.1 Aerodynamic Derivative Generation - The rigid 
body aerodynamic data for PAST I was generated from wind 
tunnel tests and Air Force DATCOM methods (reference 6) . 
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Data on the F-8 SCW aircraft (reference 2) and the basic 
Firebee II drone (reference 3) were compared to the DAST I 
estimated values where applicable to generate con- 
fidence in the estimates. No flight manuevers are planned 
to collect performance or stability and control data unless 
significant discrepancies are encountered. 

The 1/6 scale DAST I wind tunnel mode is shown in fig- 
ure 3.5. Static derivatives from M= 0.4 to M = 0.98 were 
determined from tests on this model in the NASA LaRC 8- foot 
Transonic Dynamics Tunnel. The static derivative values 
were extended to M = 1.2 by extrapolating from trends ob- 
served in the F-8 SCW flight test results. 

The static longitudinal data is presented in Appendices 
A. 1.1 through A. 1.3. Some general observations about the 
static longitudinal derivatives are summarized below: 

1. The lift, drag, and pitching moment coeffic- 
ients are non-linear in both angle of attack 
and Mach number. 

2. The value of Cl can exceed 7.7 per radian at 

a 

small angles of attack in the vicinity of 
M = 1.0. 

3. A severe pitch-up characteristic can be ex- 
pected because of a large positive value for 

C m above approximately seven degree angle of 

a 

attack. 
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FIGURE 3.5 OAST I ONE-SIXTH SCALE WIND TUNNEL MODEL 




4 . 


The control surface effectiveness in the 
pitch axis at high Mach numbers is approxi- 
mately twice the value seen on the basic drone 
configuration. 

Some general comments on the static lateral-directional 
derivatives found in Appendices A. 1.4 through A. 1.6 are: 

1. The yawing moment coefficient, C n , is non- 
linear with respect to sideslip angle. 

2. The rolling moment coefficient, C^, experiences 
rather large non-linearities at high angles of 
attack (above M = .80). 

3. Non-linearities with angle of attack exist in 
C n - 

4 . Significant non-linearities exist in both Mach 
number and angle of attack for: 

• c n g ' c £ g 

• C *«R 

• c y<s a ' Cn 6 a 

5 . Derivatives non-linear with Mach number alone 
are: 


The generation of the dynamic derivative estimates were 
accomplished using Air Force DATCOM procedures. The limit- 
ing Mach number for which the estimates are valid is M = .90. 
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The higher Mach number values were extrapolated from F8-SCW 
and Firebee II trends. The static derivative values gener- 
ated in the wind tunnel tests were used in the DATCOM equa- 
tion whenever appropriate. 

The estimated dynamic derivatives and the equations 
used to determine them are presented in Appendix A. 2. 

General comments observed in this simulated data are: 

1. Because the static derivatives were non- 
linear in either Mach number or angle of attack 
(or both) , these trends appear in the dynamic 
derivatives as well. 

2. The higher aspect ratio ARW-1 and increased 
sweep angle of the wing causes a considerable 
improvement in the longitudinal dynamic deriva- 
tives compared to the basic Firebee II. 

3. Those lateral-directional derivatives where the 
vertical tail is the dominant contributor are 
very similar to the Firebee II values as should 
be expected. 

4. The derivatives C n and Co show considerable 

P P 

improvement from the Firebee II values due to 
the increased aspect ratio. 

5. Yaw damping, C nr , is considerably decreased in 
comparison to the Firebee II , again due to the 
increased aspect ratio of ARW-1. 

6. The variation in C with angle of attack is 
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much greater for the DAST I than for the Firebee II. 

Some confidence is obtained from these estimates when 
it is observed that the non-linearities that are evident in 
the DAST data were also observed in the F8-SCW . 

3.2.2 Stability and Control Characteristics - The open 
loop flying characteristics for DAST I were defined based 
on the data presented in the previous section. The equations 
used in the longitudinal analysis were: 

[cj {i} - M |x| + [b] {u} (3.1) 

where 
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The following center of gravity corrections have been 
included in equations (3.3), (3.4), (3.7), and (3.8). The 
prime terms represent the correct values: 
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(Continued) 


(3.10) 
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The transformations listed above were exercised when- 
ever the center of gravity differed from the 25 percent MAC 
that the aerodynamics were referenced to in the wind tunnel 
tests . 

The lateral-directional analysis was done in the stab- 
ility axis system. The lateral-directional aerodynamics 
in the appendix are referenced to the body axis coordinate 
system. The axis transformation equations found in Appendix 
A of reference 7 were used to place the inertias and deriva- 
tives into the stability axes. 

The center of gravity envelope for DAST I is shown in 
figure 3.6. This plot shows the movement in the center of 
gravity with fuel burn-off. The forward and aft center of 
gravity limits are based on recommendations in reference 5. 
If the center of gravity of DAST I remains within this 
range, the loads produced by the horizontal tail will be 


35 


GROSS WEIGHT, KGS 


1085. OOi 


1040.00 


• 995.00 -I 


LC 

(S 


u: 

u. 


c c 

c 

r 

<\ 

V. 



FUSELAGE STATION 
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within the range of those produced with the basic Firebee 
II drone. The target zero weight center of gravity is re- 
commended to fall between fuselage station 269.0 and 270.0. 
This allows a margin of safety at the full fuel and zero 
fuel extremes. 

The predicted inertias and engine moment arm in pitch 
are shown in Table 3.3. These estimates are based on in- 
formation in reference 5. The engine moment arm is a meas- 
ure of distance between the vertical center of gravity and 
the engine thrust line. 

Trim conditions at various locations in the flutter 
testing region were established in order to determine the 
rigid body open loop roots. Requirements for throttle set- 
ting, elevator deflection, and trimmed angle of attack 
throughout the testing region are shown in figure 3.7. These 
estimates indicate that the angle of attack throughout the 
flutter testing region will not exceed four degrees which 
is well below the pitch-up condition. Trim elevator re- 
quirements also will remain relatively small. Based on the 
engine model used, excess thrust is still available at the 
low altitude, high Mach number extreme of h = 3048 meters 
and M = 0.98. However, the 80 percent minimum throttle set- 
ting prevents some low speed testing at lower altitudes . 
Operation of the engine below the 80 percent setting is 
avoided because of engine overheating. 
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TABLE 3.3 


ESTIMATED INERTIAS AND CG LOCATIONS 



Full Fuel 

Zero Fuel 

*xx 

259.51 

186.37 

1 yy 

3277.89 

3105.83 

Izz 

4194.62 

3266.49 

-^xz 

19.54 

16.55 

& 
o 
1 X 

13.4 

28.7 

z cg 

.074 

.146 


* Inertias are in Kgs - m 2 
cg's are in percent MAC 
MAC = .68707 meters 
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The unaugmented aircraft roots were determined by im- 
plementing the DAST rigid body equations of motion in the 
CONTROL computer program of reference 8. The short period 
mode locus is shown in figure 3.8. From figure 3.8 the worst 
and best short period characteristics can be determined. 

The worst condition occurs at an altitude of 7620 meters 

and M= 0.98. Here the short period damping ratio, c , is .15 

S P 

with a natural frequency, w n , of 8.19 rad/sec. This is at 

S P 

the full fuel condition. The best characteristics occur at 
t ' 

h=3048 meters and M=0.9 with no fuel left. At these conditions, 

^ sp = * 372 and w n =5.48 rad/sec. The short period frequency 
sp 

increases rather substantially between M=0.9 and M=0.98 due to 
large increases in C 

m 

a 

The phugoid mode locus is presented in figure 3.9. 

Here it can be seen that the worst case condition occurs 
with no fuel at M=0.8 and h=3048 meters. Under these circum- 
stances, the phugoid natural frequency is .0473 rad/sec with 
a damping ratio of -.0173. This corresponds to a time to 
double the amplitude in the phugoid mode of T 2 =847 seconds. 

P 

The best phugoid mode characteristics occur at an altitude of 

3048 meters and M=0.98, again with no fuel on board. At these 

conditions w =.0302 rad/sec and C =.0484. The changes in the 
X P P 

phugoid mode characteristics with increasing Mach number are 

due to a change in sign of c from negative to positive. 

u 

In the lateral-directional case, the dutch roll mode 
locus is plotted in figure 3.10. Minimum damping in this 
mode occurs at the conditions of full fuel, h=3048 meters. 
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and M= 0.9. At this flight condition, co =5.85 rad/sec and 

n D 

C d = 0.181. The best dutch roll mode characteristics are ex- • 

hibited at zero fuel, h=7620 meters, and M=0.7. The natural 

frequency is 4.576 rad/sec with ? D = .332. 

The roll and spiral mode loci are shown in figure 3.11. 

For the roll mode, the worst case condition is at full fuel, 

h « 7620 meters, and M = 0.7. The roll mode time constant, 

T_ , is .741 seconds at this condition. When the aircraft is 

at the zero fuel condition at M=0.98 and h = 3048 meters, then 

T= .106 seconds for the best characteristics. 

R 

With the DAST I in the empty fuel configuration, there 
are flight conditions where the spiral mode is stable. The 
most stability is achieved at M = 0.8 at h = 3048 meters. The 
spiral mode time constant, T g , is 357.14 seconds at this flight 
condition. The spiral mode is unstable at all flight conditions 
with full fuel. When M = 0.98 and h = 7620 meters, the 
worst case occurs with T = -84.03 seconds with a time to 

y 

double the amplitude of 58.25 seconds. 

To place the unaugmented characteristics of the DAST I 
in their proper perspective, the aircraft was defined as a 
Class II (reference 9) vehicle which will perform Category 
B flight phase type manuevers. This corresponds to trans- 
port aircraft requirements of gradual manuevers with accurate 
flight path control. 

Under this criterion the DAST I possesses the following 
unaugmented flying qualities according to reference 9. 
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CHAPTER 4 


CONTROL SYSTEMS SYNTHESIS 

This chapter will describe the augmentation control 
systems developed for the DAST I. A description of Boeing's 
analysis for the FSS design is presented first. The develop- 
ment of the primary RAV flight control system is presented 
second. The backup autopilot system is described third with 
some special requirements for launch and recovery examined. 

4.1 Flutter Suppression System - In the scheme called 
flutter suppression, an augmentation control system gener- 
ates artifical stiffness in a system, such as ARW-1, by 
producing forces which oppose motion in the system. Diver- 
gent motion of the wing due to flutter is designed to be 
suppressed by the FSS as long as that control system operates 
properly. 

The design of the FSS was conducted by the Boeing Wich- 
ita Company under contracts NAS1-14675 and 14028 from NASA 
LaRC. The major contract requirement was to provide a flut- 
ter suppression system to extend the flutter boundary of 
DAST I by twenty percent. Additional detail on the informa- 
tion presented in the next section can be found in reference 
10 . 

4.1.1 Equations of Motion - Through the design cycle, 
vibration solutions for the structural model described in 
section 3.1.1 were obtained. For each solution, structural 
reference points were added at the midpoint of the front 
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and rear spars at each rib and on the hinge line of the con- 
trol surface. The points, also called co-location points, 
contain the motion for the six degrees of freedom and are 
used to obtain the motion of the aerodynamic panels in the 
aerodynamic analysis. 

The flutter analysis sequence is illustrated in figure 
4.1. The unsteady aerodynamics forces on the wing and hor- 
izontal and vertical tail were generated using a bwc develop- 
ed three-dimensional plate doublet finite element solution. 
The theory accounts for Mach number and finite span effects, 
and includes the aerodynamic coupling between all the com- 
ponents of the aircraft. The finite element wing aerodynam- 
ic model is illustrated in figure 4.2. The unknown pressure 
distribution for each aircraft mode is determined by assum- 
ing the pressure to be constant over a given aerodynamic 
panel and solving for the pressure based on a specified Mach 
number and reduced frequency. The primary lifting surfaces 
are modeled with trapezoidal elements arranged in strips 
parallel to the free stream. 

The equations of motion for the flutter analysis were 
initially formed using complex oscillatory aerodynamic co- 
efficients determined for specific values of reduced fre- 
quency, uj/U 0 . The final form of the equations of motion 
were formulated in terms of real matrices through introduc- 
tion of an approximating function. 
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FIGURE 4.2 DAST ARW-1 AERODYNAMIC MODEL 


From reference 10 the original equations in standard 




generalized modal coordinate 
aerodynamic influence coefficient matrix 
the linearized boundary conditions 
free stream velocity 
vertical gust velocity 


Only the symmetrical equations will be examined in this 
report. The anti-symmetrical equations will be of the same 
form and can be assertained from reference 10 . 

The approximating functions were developed by plotting 
particular elements of the complex matrix A x at various 
values of u from 0 to 400 radians/second. This is illustrat- 
ed in figure 4.3. The solid line in figure 4.3 is an approx- 
imating function, chosen as a rational polynomial function 
of the complex variable s. The use of the approximating 
function allows accurate approximation of the time delays 
inherent in the unsteady aerodynamics subject to the fol- 


lowing restrictions: 
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• It must have complex conjugate symmetry. 

• It must have denominator roots in the left half 
plane . 

• It must approximate the value of the complex 
coefficient when s = o + jtu, for those values 
of to analyzed. 

The approximating function for each element in the Aj 
matrix was determined after analysis at seven discrete fre- 
quencies. By substituting the approximating functions for 
the aerodynamic influence coefficients, the final form of 
the equations of motion results. 


(s 2 [m] + s [Z] + [k]) <jq(s)} 
+ (s 2 p jbj + pU 0 s jc 2 j + pUj 

4 

{q(s)j+ (pU 0 [R 0 ] + fUo 



(4.2) 



The first line in equation (4.2) represents structural 
coefficients. The second line is aerodynamic coefficients 
and the third line is gust velocity coefficients. 

With this equation, the air density, p, and velocity, 

U Q , can be varied so solutions can be found at all the flight 
conditions of interest. Solutions of equation (4.2) with- 
out control system augmentation yields the results in 
section 3.1.2. 


4.1.2 Design Synethesis Techniques - The design criteria 
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established for the DAST ARW-1 FSS are: 

• Extend flutter boundary by 20 percent. 

• Design gust: 1.83 m/sec rms 

• System stability at design speed (MIL-F-9490D) 

- Elastic Modes: GM = +6.0 DB; PM = +45 deg. 

- Maneuver Modes: GM = +4.5 DB; PM = +30 deg. 

• No FSS redundancy 

• No degration in structural mode damping (other 
than the flutter mode) below Z = .01 except when 
the unagumented mode damping is less than Z = .01. 

• A single wing flutter suppression system will be 
analyzed but not implemented. 

The synthesis started with a parametric evaluation to 
determine combinations of control surface locations and sen- 
sor type and location that would control the flutter mode 
and minimize the coupling with other wing structural modes. 
This parametric evaluation was achieved by using a zero 
locus to estimate modal stability and coupling as a func- 
tion of sensor and control surface placement. If the sen- 
sor or control surface placed the zeros in the right half 
of the s-plane near the flutter mode it was eliminated. This 
would indicate the coupling was small. Other combxnatxons 
were eliminated because their location in the right half 
of the s-plane indicated strong adverse coupling with other 
structural modes. Using this procedure, along with mode 
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shape information, candidate sensor and control surface loca- 
tions were determined. Control laws were then synthesized 
and were evaluated using root locus methods. 

The control surface size and location: which was selected 
has the following characteristics: 

• 23 percent chord hingeline 

• .254 meter span 

• Inboard edge is 1.71 meters (WBL 67.35) from the 
centerline. 

• Outboard edge is 1.964 meters (WBL 77.35) from 
the centerline. 

The sensor location is 2.007 meters (WBL 79) from the 
fuselage centerline and is attached to the rear spar. 

This combination places the zeroes of the higher fre- 
quency modes very near their corresponding poles. The flut- 
ter mode zero is in the left half of the s-plane and influ- 
ences the unstable flutter mode root into the left half of 
the s-plane when the control law in figure 4.4 is mechanized. 
This is the control law for the DAST FSS as of March, 1978. 
Some adjustments from the system shown in figure 4.4 are 
possible prior to delivery to NASA. 

From figure 4.4 it can be seen that four accelero- 
meters are utilized. One on each wing at the location spec- 
ified before, one on the fuselage centerline at FS 265 and 
a roll angular accelerometer on the fuselage centerline at 
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FS 274. The last two accelerometers on the fuselage center- 
line are used to subtract the rigid body vertical and roll 
accelerations from the wing sensors. The filter after the 
second summing junction has a parameter, D, in numerator and 
denominator terms which is gain scheduled with Mach number 
and altitude. Its purpose is to maintain adequate gain mar- 
gins at both low and high dynamic pressure conditions. A 
washout filter (s/s + 2) eliminates the influence of the 
very low frequency rigid body modes. The numerator term, 

(s + 370), is compensation for an actuator pole. The sec- 
ond order lead term (s 2 + 300s + 600 2 ) is to achieve the 
proper high frequency phase and gain characteristic, as 
are the lag terms (s + 40) 2 and (s + 800) . The lag terms 
(s + 1000) and (s + 2000) are for high frequency attenua- 
tion. Notch filters for both symmetrical and anti- symmetric- 
al conditions are also seen in figure 4.4 and are to achieve 
the proper gain and phase at the corresponding flutter fre- 
quencies. 

A root locus of this control law at the conditions of 
M = .9 and 3048 meters is shown in figure 4.5. Boeing's 
analysis showed that about 90 degrees of phase lag was re- 
quired to damp the flutter mode. This amount of phase lag 
is typical of using an accelerometer for the feedback sen- 
sor giving a pseudo rate signal at the flutter mode fre- 
quency. It also provides attenuation of the higher f re- 
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quency modes thus helping to decouple them from the FSS. 

The surface deflection limit for the wing FSS aileron 
is + 12 degrees. The aileron surface displacement and rate 
requirements are shown in figure 4.6. The requirements for 
3.67 m/sec peak gusts are a surface deflection of 18.135 de- 
grees and a surface rate of 1460 deg/sec. Clearly with a 
12— degree deflection limit, saturation of the control sur- 
face will occur. Increasing the surface size to eliminate 
this problem was not an option since increased hydraulic 
power requirements could not be met with the pump restric- 
tions that existed for this vehicle. 

4.1.3 Augmented Flutter Boundary - The improvement 
in the open loop flutter boundary which the FSS makes is 
shown in figure 4.7. If the jettisonable mass ballast at 
the trailing edge near the tips is removed, then the new 
open loop flutter boundary is shown. As long as the DAST 
vehicle is below this second boundary and encounters dif- 
ficulty with the FSS, jettisoning the ballast will provide 
some improvement in safety. However, if the vehicle is 
flown above this boundary and difficulty is encountered, 
then even if the mass balance is jettisoned, flutter will 
occur. 

To avoid exceeding M = .98 and still achieve the twenty 
percent increase in the flutter boundary , the DAST vehicle 
must fly at low altitudes. An altitude minimum for safety 
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reasons has been set at 3048 meters. The improvement in 
the flutter mode damping ratio with the FSS operating is 
illustrated in figures 4.8 and 4.9. It can be seen that the 
predicted inprovement in flutter mode damping exceeds the 
required twenty percent goal. The BWC plate-doublet aero- 
dynamics were used through M = 1.1 in the analysis. The 
accuracy of a linear subsonic aerodynamic theory to this 
high a Mach number is questionable but BWC maintains that 
their prior experience in this area shows that although the 
theoretical results will be less accurate than at lower Mach 
number, they will be conservative. Thus their use in the 
analysis is considered appropriate. 

4.2 RAV Primary Control System . - The remotely augmented 
control system will be the primary flight control system for 
the DAST RPRV. This system which was illustrated in figure 
2.3 allows great flexibility in control system options. The 
DAST I will exercise a three mode RAV system. By this, it is 
meant that the pilot has the option of selecting three differ- 
ent control system configurations at his discretion. The 
three modes available for DAST are: 

1. Computer Direct (CD) 

2. Rate Damper (RD) 

3. Command Augmentation (CAS) 

The computer direct mode does utilize the RAV ground 
computer to determine pilot stick commands but does not aug- 
ment the basic aircraft. The rate damper mode adds rate 
dampers to the pilot inputs. The damper gains are pilot 
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selectable. The command augmentation mode modifies the pilot 
inputs dependent upon flight conditions. The stick commands 
are input into an integrator, whose characteristics can be 
changed by certain flight conditions. In the CAS mode, the 
pilot does not have direct control of the control surfaces. 

4.2.1 Design Considerations . - The RAV control system 
requirements were: 

• Augment the DAST I to level 1 flying qualities 

(MIL-8785B) in all axes at all flight conditions. 

• Control the unstable pitch-up characteristics 
of the unaugmented vehicle. 

• prevent the vehicle from exceeding the 2.5 g 
wing load limit. 

The above requirements were established in order to reduce 
the pilot's workload during normal and emergency conditions. 
By offering load protection and controlling the vehicle's 
abrupt pitch-up tendency above approximately seven degrees 
angle of attack, the pilot can perform maximum performance 
pull-up maneuvers to reduce speed as fast as possible should 
an emergency arise in supercritical FSS testing. 

The RPRV philosophy of no redundant systems was a factor 
in determining the type of augmentation supplied by the RAV 
ground computer. Recall that in section 2.2.2, the Air 
Force requirement for a secondary command capability was 
established. However, this backup system and the RAV system 
share many common and non-redundant systems . This considera- 
tion led to the establishment of the following design guide- 
lines : 
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• Systems common to the RAV and on-board backup auto- 
pilot systems (servoes, hydraulic power supply, 
electrical power supply, etc.) have significantly 
better reliability than the mid-air recovery (MARS) 
technique, thus a single failure of common systems 
may result in the loss of the vehicle, but its oc- 
curence is much less likely than a MARS failure. 

• The RAV and on-board control systems will have no 
sensors in common from the single sensor set in order 
to minimize the number of common systems. 

Using these guidelines guarantees that no sensor failure 
can cause a degradiation in the performance of both control 
sys terns . 

Due to the role which the autopilot performs, the attitude 
gyros and angle of attack sensor were reserved for its use. 

This is covered in detail in section 4.3. With the restriction 
from use of these sensors in the RAV system, the synthesis of 
the primary flight controls evolved into the system described 

in sections 4.2.2 and 4*2.3. 

4.2.2 Longitudinal Modes . - A block diagram of the RAV 
system implemented in the pitch axis is shown in figure 4.10. 

As mentioned previously, the CAS mode will be the primary op- 
erating mode for the flight tests. It should be noted that 
the flight test experiment could be accomplished in the RD 
mode but the pitch-up and load limiting protection offered in 
the CAS mode would not exist. To fly the vehicle under these 
circumstances increases the risk and the pilot workload in- 
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creases substantially. The philosophy behind the multi-mode 
option is twofold. First, if the pilot encounters flying 
characteristics significantly different from those predicted, 
the RD mode will be less sensitive to these changes and will 
give the pilot direct control of the vehicle's control sur- 
faces, instead of operating through integrators as in the CAS 
mode. Second, the failure detection software described in 
section 5.5.4 can detect dangerous sensors failures and down- 
mode to the next lowest mode which does not require that sen- 
sor . 

The limiter and fader block in figure 4.10 serves the pur- 
pose of putting the integrator in a hold status above the 
2.5 g load limit. A trigger g level of 2.35 g's first acti- 
vates the limiter and is required due to the fact that the 
aircraft’s response lag's the pilot's inputs. The action of 
placing the integrator in hold prevent's the pilot from con- 
tinuing the pull-up and over stressing the wing. The fader 
simply reduces the gain on the integrator between the 2.35 
g trigger level and 2.5 g's so that continued pilot back 
pressure on the center stick allows the maximum 2.5 g man- 
euver without overshoot. 

Normal acceleration feedback above 2.5 g's provides pro- 
tection in turbulence when large g maneuvers are performed. 

A lead- lag filter on the pitch rate signal was required 
to compensate for the one sample delay in the control law 
computation. The one sample delay is caused by the finite 
time interval required to make the control law computation. 
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This requires that past information must be used in the most 
current computations. Thus, introducing a time lag which is 
dependent on the age of the data used in the control laws. 
Without the lead, the additional delay in the pitch rate at 
high Mach number results in a signal similar to angle of 
attack (lagged pitch rate) . Feeding this back only increases 
the short period frequency without improving the damping 
ratio. 

Since the RAV system is a sampled data system, z-plane 
root locus methods were used in its design. The CONTROL 
program of reference 8 was used for the analysis. The fly- 
ing qualities are defined in terms of the s— plane character- 
istics so the z-plane information can be related to the s- 
plane by the relationships: 


U). 


arctan y 


(4.3) 


5 = COS 


arctan 


-arctan y 


‘■n 


cos arctan ~ 


(4.4) 


where v is the imaginary part of the z-plane complex root, 
y is the real part, and T is the sample period. 


A z-plane root locus of the longitudinal RAV system is 
shown in figure 4.11. Lines of constant damping ratio appear 
as logarithmic spirals in the z-plane. These spirals are 
shown in figure 4.11. The flight condition corresponds to 
the worst case short period characteristics. Open loop C gp 
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= .15, while in the RD mode, is increased to . 386. In 

sp 

the CAS mode the short period damping ratio is .356. This 
occurs at a feedback gain of .3 deg/deg/sec. The pitch CAS 
integrator is a pure integrator and therefore, its gain re- 
mains 1.0. The phugoid damping ratio, although not the 
worst case, is representative of how the CAS system augments 
it. Its open loop value in this case is 0.015, in RD s pH =.016, 
and in the CAS mode with K = .3, c pH is greater than one. 

Under the best circumstances, that is h = 3048 meters and 
M = .90, and with the same gains, c g p * .774. Since only trans- 
port-type maneuvers will be performed with DAST I, a washout 
filter in the pitch rate feedback path was deemed unnecessary. 
The phugoid damping ratio at this condition is greater than one 
in the CAS mode. 

To show that the pitch RAV system is also adequate for 
unplanned emergency conditions, a pull-up maneuver was examined 
which is at the approximate condition where the pitch-up char- 
acteristic is evident. This might occur as the vehicle is in 
a maximum performance slow down and climb out following an 
FSS problem. As the vehicle slows to M = .7 going through 

h = 6096 meters, C begins to become positive since the 

a 

vehicle is at approximately seven degrees angle of attack. 

For the open loop aircraft, its short period roots have migrat- 
ed to the real axis, once the angle of attack approaches seven 

degrees. When C becomes positive, the short period mode 

a 

has two real roots with one outside the unit circle. This 
situation is remedied with the pitch CAS system as can be 
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seen in figure 4.12. The short period frequency becomes 2.54 
rad/sec with a damping ratio t = .537. The z-plane locus in 

^Xr 

in figure 4.12 is of the normal acceleration feedback gain. The 
gain for normal acceleration feedback was found to have a sig- 
nificant effect on the short period frequency and damping . A 
gain of .2 deg/g proved to be the most satisfactory from sim- 
ulation studies. 

In order to compensate for any aeroelastic effects or 
other uncertainties the vehicle might exhibit in flight, a 
sensitivity analysis on the dominate longitudinal derivatives 
was conducted. A + 20 percent variation in each derivative 
was examined for any adverse impact. For the short period 
mode the following effects were observed: 


Changes in C L , C m , and C m had little affect on the 
a q a 

short period frequency 


Increasing C by 20 percent increased <o n 

a sp 


at the 

worst case conditions by approximately 15 percent 
Short period damping was affected in the following 
manner at the worst case conditions: 

• is unchanged with a +20 percent change in 

C m. 

a 

• decreases by 5*6 percent with a -20 percent 
change in C L 

a 

• £ decreases by 6 percent with a +20 percent 

sp 


change in C 


m 

a 

£ decreases by 2.8 percent with a -20 percent 


change in C 


m 
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FIGURE 4.12 NORMAL ACCELERATION Z-PLANE ROOT LOCUS 
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With these worst case derivative changes added to the 
rigid body conditions existing at the worst case flight con- 
ditions (which is very unlikely to be encountered) , the aug- 
mented short period characteristics become oo = 10.55 rad/ 

sp 

sec and ; gp = .305. These are still within level 1 flying 
quality ratings. The pitch CAS will satisfy its requirements 
for the short period mode at all flight conditions in the 
flutter testing region without the requirement for scheduling 
feedback gains with dynamic pressure. 

In the case of the phugoid characteristics , which are not 
as critical, the speed derivatives were varied by the +20 
percent margin. The effects are summarized below at the worst 
case conditions for the phugoid mode: 

• Increasing C by 20 percent decreased o by 

m u n PH 

six percent and s pH by four percent 

• Increasing C- did not change o but did decrease 

T u PH 

5 pH by 2.5 percent 

• Decreasing C n did not effect w but did decrease 

D u PH 

C_.„ by two percent 

rn 

• Decreasing C T decreased w by .8 percent but did 

L U n PH 

not effect c pH 

Adding the worst case conditions listed above to the aug- 
mented values still results in a phugoid damping ratio orders 
of magnitude larger than MIL8785B Level 1 requirements. 

4.2.3 Lateral-Directional Modes - The predicted open-loop 
lateral-directional characteristics for DAST I indicated that 
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the lateral-directional RAV system is less critical than the 
pitch axis. This is because in the flutter testing region, 
the open loop aircraft exhibits Level 1 flying qualities al- 
ready. The dutch roll mode for the DAST is predominated by 
a rolling oscillation. Piloted simulation showed that because 
the roll inertia for this vehicle is so small, pilots preferred 
more dutch roll mode damping than the aircraft exhibited. 

Roll response was considered satisfactory. 

The improvement in the dutch roll damping was achieved with 
the control system shown in figure 4.13. No augmentation was 
applied in the roll axis. Since the roll response was already 
adequate, a roll damper was not necessary nor desirable 
since it had a destabilizing effect on dutch roll 
damping . 

The yaw axis utilizes washed-out yaw rate feedback in the 
RD mode. It was not necessary to add lead-lag to the yaw rate 
signal because the dutch roll mode frequency was not high enough 
to be affected by the .01875 second sample delay. Since the 
airplane has a relatively large diherdral effect, the yaw 
damper was found to be very effective in augmenting the dutch 
roll mode damping. 

2 — plane root locus of the lateral— directional RAV control 
system is shown in figure 4.14. This flight condition cor- 
responds to the worst case open loop characteristics of the 
dutch roll mode. The yaw damper increases the dutch roll 
natural frequency from 5.85 rad/sec to 6.89 rad/sec and aug- 
ments the damping ratio from ? D = .181 to S D = .449 at a 
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FIGURE 4.14 LATERAL-DIRECTIONAL RAV Z-PLANE ROOT LOCUS 


gain of .3 deg/deg/sec. Under the best open loop circum- 
stances, the yaw damper achieves a damping ratio of 0.596. 

The lateral-directional RAV system increases the roll time 
constant very slightly and the spiral is unaffected. 

Due to the large amount of damping offered by the yaw damp- 
er, the sensitivity analysis of the dominate stability deriva- 
tives in the lateral directional equation was not of concern. 
The worst circumstances for the dutch roll mode damping show- 
ed a 13 percent reduction in c D . At the worst flight condi- 
tions for the roll mode, the time constant is increased by 
approximately ten percent. Only the derivative C n seemed to 
have an appreciable effect on the spiral mode. It increased 
the spiral mode time constant by only 6.5 percent. It should 
be recognized that the reason the +20 percent variation in the 
lateral-directional stability derivatives was not of great 
concern was that, even with the worst conditions , the RAV 
system implemented for the DAST I possesses damping and re- 
sponse characteristics substantially better than the MIL- 
8785B Level 1 minimums in the lateral-directional axes. 

4.2.4 Control Law Functional Validation . - The RAV control 
laws were implemented in the DAST I real-time simulation. The 
z-plane transfer functions were coded in FORTRAN as difference 
equations. The RAV performance was tested in all areas of 
the flight envelope with a .3048 meter/second RMS gust in all 
axes . 

Figure 4.15 shows time histories of the longitudinal para- 
meters for the open loop (CD) and augmanted (CAS) at normal 
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one g flight and at the maximum 2.5 g condition. A smooth- 
ing of the angle of attack, pitch rate, and normal accelera- 
tion traces is evident when the augmented system is activated 
in the one g condition. Traces without the gust inputs verified 
the damping and frequency values resulting from the root locus 
analysis. Figure 4.15 includes the turbulence to illustrate 
the effect the RAV system will have in a more realistic en- 
vironment. The high g condition, also shown in figure 4.15, 
represents a maximum performance pull-up maneuver which might 
be initiated if problems occur with the FSS. It basically 
illustrates the proper operation of the load limiting cap- 
ability of the pitch CAS. The pilot pulls back on the stick 
and holds it against the 2.5 g limit. As his back pressure 
on the stick increases, the fader allows the pitch orientation 
integrator to slowly increase the surface deflection since 
the change in flight condition requires more deflection to 
maintain the maximum 2.5 g value. The fader's purpose, which 
is demonstrated here, is to reduce the gain on the integrator 
to a value (1/20 of normal) which keeps the vehicle at its 
maximum performance load limit without overshoot. The nor- 
mal acceleration feedback above the 2.5 g limit is used to 
protect the vehicle from gust inputs at the limiter value. 

Figure 4.16 shows the operation of the RAV lateral-direc- 
tional system at the same one g flight condition used in 
the longitudinal case. The rudder and tail doublets during 
operation in the computer direct mode (CD) illustrates the 
following: 
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FIGURE 4.16 TIME HISTORY OF THE LATERAL-DIRECTIONAL 
RAV OPERATION 
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• The dutch roll mode involves much more rolling 
oscillation than yawing motion. 

• The rolling tail has little effect in exciting the 
dutch mode. 

• Rudder doublets are very effective in exciting the 
dutch roll mode. 

This confirms the reasoning for augmenting the dutch roll 
mode damping with the yaw damper and not using a roll damper. 

With the yaw damper operating, a significant reduction in 
the amplitude of the roll rate is evident. Doublets to the 
rolling tail and rudder illustrate the improvement in the 
dutch roll mode achieved with the lateral “directional RAV 
system in a gusty environment. 

Pilot evaluations of this vehicle drew the following com- 
ments : 

1. The pilots were surprised by the amount of improve- 
ment in flying qualities resulting from the ARW-1 

wing . 

2. All the control modes perform their desired tasks 
satisfactorily. 

3. The pitch CAS mode reduces the pilot workload parti- 
cularly in turns. 

The third comment has significance in that the pilot felt 
that the pitch CAS mode allowed for time to concentrate on 
maintaining a constant airspeed which is critical in the 
flutter tests. 

4.3 PAST I On-Board Autopilot. - The basic Firebee drone 
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had an on-board discrete command anlaog autopilot system. In 
order to interface this system properly with the digital 
proportional RAV system, a complicated interface box would 
have been required. For the requirement as an RPRV, only 
approximately 10 percent of the Firebee autopilot's function- 
al capabilities were utilized. The other autopilot capabili- 
ties were associated with its role as a military target drone. 

A new much simpler on-board autopilot system was developed 
for the following reasons: 

1. The cumbersome Firebee autopilot and its interface 
equipment required a large portion of the room 
available in the DAST fuselage. 

2. It was desirable to simplify the mechanization of 
the control systems by making the autopilot propor- 
tional like the RAV system. 

3. Modifications to the Firebee autopilot to adapt it 
to the DAST flying characteristics would have been 
difficult since the Firebee autopilot construction 
was of multi-layer printed circuit boards. 

4.3.1 Design Requirements . - The primary purpose of the 
DAST I autopilot is to take over command of the vehicle if 
the RAV control system malfunctions. It will be used only to 
fly the DAST safely to the MARS recovery area. With this in 
mind, the following requirements were established: 

• Before F-104 chase captures command: 

-wing leveling 
-envelope limiting 
* angle of attack containment 
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* limit velocity by a-control 
-Have at least Level 3 flying qualities (with 
+20 percent uncertainty in aerodynamic derivatives) . 

• After F-104 chase captures command: 

-bank angle control 

-glide slope control with at lease Level 3 
minimums 

-Throttle control 

If there is a loss of TM which causes RAV failure , the 
autopilot is to stabilize the vehicle at the last heading and 
altitude until the F-104 takes control or an automatic recovery 
is initiated. The throttle setting is automatically reduced 
to the 80 percent RPM setting when there is an RAV-to-AP 
transfer. Since the chase aircraft RPRV pilot cannot tell 
the magnitude of angle of attack at which the DAST is flying , 
the angle of attack must be limited to prevent inadvertant 
high angle of attack operation and pitch-up. 

Priority in sensor selection was given to the autopilot, 
since the RAV computer gave the ground facility the capability 
to synthesize any parameter which it required. 

4.3.2 Longitudinal Axis . - A block diagram of the complete 
autopilot is shown in figure 4.17. Pitch and roll commands 
are sent to the vehicle in the form of left and right horizon- 
tal stabilizer commands. These are symmetric for elevator 
inputs and differential for roll inputs. In this manner, 
the same digital word format used in the RAV uplink is used 
in the AP uplink. These digital commands are converted to 


92 




FIGURE 4.17 DAST I AUTOPILOT BLOCK DIAGRAM 
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analog signals on-board since the autopilot is an analog 
system. 

The longitudinal functions in this system are: 

• pitch damping 

• pitch attitude hold 

• angle of attack limiter 

• glide slope hold 

Pitch damping was achieved by differentiating the pitch 
attitude signal to produce an approximate pitch rate signal. 
Pitch attitude feedback augments the phugoid characteristics . 

An s-plane root locus of the pitch damper and pitch at- 
titude feedback loops is presented in figure 4.18. The 
flight condition corresponds to the worst case short period 
open loop characteristics in the flutter testing region. The 
short period frequency is increased to 9.56 rad/sec and the 
damping ratio increased to .208. The phugoid frequency 
jumps to 0.044 rad/sec while the damping ratio is increased 
from .015 to .79. The feedback loop gains are .05 deg/deg/ 
sec for the damper and .05 deg/deg for the pitch attitude 
hold loop. The short period gain was chosen to provide suf- 
ficient margin so that when the +20 percent variation in the 
aerodynamic derivative values are included, at least Level 3 
minimums exist. The high phugoid damping value was found to 
be desirable from simulation studies. 

By controlling the vehicle's glide slope, both velocity 
and altitude excursions can be kept within safe limits. For 
the DAST autopilot, the glide slope signal, y, was generated 
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by subtracting angle of attack from pitch attitude. Since this 
signal is not bank angle compensated, the value for y will 
only be accurate when the bank angle is not large. A root 
locus of the glide slope loop is shown in figure 4.19. The 
forward loop integrator is required to magnify, over a period 
of time, the error between the glide slope feedback path and 
the reference command so the steady state commanded y can be 
achieved. The lag filter for the angle of attack signal was 
required to eliminate the nose boom structural dynamics from 
the signal. Due to the low frequency nature of flight path 
control, the glide slope hold has little effect on the ve- 
hicle's short period mode. With a pure integrator and a y- 
loop gain of .1 deg/deg, the phugoid damping ratio becomes 
.767 at a frequency of .112 rad/sec. 

At Mach numbers where the phugoid mode is unstable, the 
pole at the origin is attracted toward a zero in the right half 
of the s-plane near the origin. At the gains mentioned 
above , a pole-zero cancellation takes place in the y/$ e 
transfer function to effectively eliminate the divergence. 

The zero which is in the right half of the s-plane can be 
approximated by: 


1 s _W_ f iH _ 3T*| (4.5) 

T Y1 g L 9U 9U J 

This equation is given to illustrate the source of this 
zero. Since the pitch attitude is tightly controlled by the 
inner loops, no divergence is evident in the pitch response. 
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This zero does not exist in the u/6 e transfer function, so 

no pole- zero cancellation occurs. The divergence shows up 

as a speed stability problem created by the fact that C D is 

negative and C T is positive in many regions of the flight 
u 

testing envelope. It makes precise control of airspeed very 
difficult. In general, the conditions which possessed an 
unstable open loop phugoid mode also experienced this speed 
stability annoyance. Since speed brakes do not exist and 
throttle feedback is impractical due to very poor engine re- 
sponse characteristics, the 1/T zero cannot be moved to 

Y i 

solve this problem. 

The alpha limiter in the longitudinal autopilot only works 
above 6 degrees angle of attack and below -4 degrees angle 
of attack. Between these limits no feedback occurs. The 
case which was set up to examine the limiter was: 

M - .8 

h = 8680 meters 
load factor = 2.5 g's 
“trim “ 7,25 degrees 

This might correspond to conditions in an emergency pull 

up maneuver to slow the vehicle down. With M = .8 at 7.25 

degrees angle of attack, C m = .241/rad which would represent 

a 

a severe pitch- up tendency. The root locus of the alpha lim- 
iter is seen in figure 4.20. Because the glide slope hold is 
outer loop in nature, the forward loop integrator had to be 
disabled so that it would not overpower the angle of attack 
limits. With the rest of the autopilot augmentation operating, 
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the short period mode has a frequency of .775 rad/sec and 

a damping ratio of t = -.127. The phugoid has become two 

real negative roots at this condition. With a gain of 1 deg/ 

deg on the angle of attack feedback, u> n = 5.02 rad/sec and 

r = .08. Lower gains on the alpha limiter increased the 
sp 

damping ratio, t, , but simulator studies showed that lower 

sp 

gains were inadequate in dynamic tests. An authority limit 
on the alpha limiter of 1.5 degrees of elevator command gave 
it the characteristics of a soft limit which eliminated the 
limit cycle tendency this high gain produced across the alpha 
boundary. 

4.3.3 Lateral-Directional Axes . - The inherent lateral- 
directional characteristics have been shown to be quite good 
for the DAST I vehicle. The autopilot design requirements 
specified only the need of a wing leveling capability. Upon 
examining a single wing leveler system using bank angle for 
feedback, it was found that no additional augmentation would 
be required to satisfy the design specifications. A root 
locus evaluation at the worst case conditions for the dutch 
roll mode is shown in figure 4.21. This mechanism increases 
the dutch roll frequency by only .16 rad/sec and increases 

r from .181 to .21. The roll mo'de time constant is increased 
by 15 percent. The spiral mode is stabilized with a time 
constant of 5.15 seconds. 

4 . 3.4 Autopilot Control Law Functional Validation . - The 
F-104 RPRV control of the DAST autopilot in the pitch axis is 
demonstrated in figure 4.22. A climb maneuver followed by 
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a straight and level command is shown. The turbulence used 
in the RAV time histories is included as a disturbance. It is 
obvious from the trace of elevator deflection that only very 
gradual maneuvering is possible in the AP mode. If no inputs 
are made, the vehicle will maintain very precise control of 
altitude. 

The function of the wing leveler is illustrated in figure 
4.23. The roll rate trace shows very little difference 
from the open loop aircraft response indicating the dutch 
roll mode is influenced very little. The turn maneuver shows 
that the bank angle is held relatively constant in a turn 
and with wings level. The bank angle command was limited be- 
cause the glide slope signal does not include the cos<j> term 
in the expression a *= (0 — a)cos$. Without this compensation, 
large errors are created in the a values when <j> becomes large. 
This error can drive the vehicle to the alpha limit creating 
high g's in the process. 

Pilot evaluations of the autopilot in the real-time simu- 
lator produced the following comments: 

1. The autopilot performs its task satisfactorily. 

2. With the present trim rates on the F-104 DAST control 

stick, adequate turn capability exists and adequate 
climb and dive performance exists. 

3. The slow response does provide the necessary load 
limiting protection. 

4 . 3.5 Launch and Recovery Requirements . - At present, the 
launch requirements for the DAST I are to launch with the 
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control surface locked in a nose down command for three sec- 
onds before the control systems become active. To examine 
the reasonability of this requirement, the flowfield effects 
of the B-52 launch vehicle were added to the real time sim- 
ulation. The flowfield model was developed for the 3/8 
scale P-15 RPRV of reference 4. This same flowfield model 
was used in the basic Firebee simulation and compared reason- 
ably well with flight results. 

The DAST I vehicle is attached to the B-52 pylon with a 
-5 degrees angle of attack relative to the B-52 wing. The 
flowfield is estimated to add another -1.8 degrees to this 
value. The vehicle will sense a -2.7 degrees sideslip angle 
when attached to the pylon. Figure 4.24 shows a simulated 
launch from 4572 meters at M ■ .65. The B-52 wing angle-of- 
attack was estimated to be 5.8 degrees. The vehicle rolled 
to the right approximately 40 degrees before the RAV system 
was activated. No unsafe conditions were encountered. 

In figure 4.25, a launch in the AP mode is shown. This 
would result if the DAST was launched active or the RAV fail- 
ed to operate after launch. About 1.5 g's of normal accelera 
tion is reached in the pull-up maneuver to achieve zero glide 
slope. The vehicle still rolls to about 40 degrees of bank 
angle and very little altitude is lost. 

The recovery technique was investigated to see if the low 
g wing would encounter problems with the MARS recovery . Esti 
mates for the drouge and main parachute drag characteristics 
were included in the equations of motion. The proper chute 
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attachment points were included so some guess at the pitch- 
ing moments caused by the parachute drag could be made. The 
drag created from the chutes was ramped into the equations 
to represent the time required to "unfurl" the chutes. The 
Firebee recovery logic was programmed into the simulation in 
the same fashion in which it interfaces on-board the vehicle. 
The vehicle has to be in the AP to recover. Time histories 
of the longitudinal parameters during a normal recovery ini- 
tiated at M = .93 and h = 3658 meters are seen in figure 4.26. 
The aerodynamic data tables only go to 14 degrees angle of 
attack so the parameters are invalid above this value. The 
dynamic pressure is so low at this point that the main chute 
deployment creates no appreciable forces as the vehicle rotates 
to 90 degrees about the main chute attachment on the wing 
center section. 

Two things were evident in the recovery study. First, if 
an emergency recovery is commanded at these same conditions, 
the wing would be overstressed when the main chute deployed. 

In an emergency recovery the main chute deploys 2.5 seconds 
after the drag chute which comes out immediately. This gives 
insufficient time to slow down before the vehicle rotates at 
main chute deployment. The resulting high dynamic pressures 
and high angles of attack would cause damage to the wing. 

The second item resulted in a modification to the basic 
Firebee recovery logic. It was found that the DAST I would 
be hard against the alpha limiter with little control surface 
effectiveness below M * .4. The basic drone would not deploy 
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TIME, SEC 

FIGURE 4.26 RECOVERY TIME HISTORY 
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the drag chute until it slowed to M s .3 in the normal re- 
covery climb out ( h> 4572 meters) . Changing this Mach num- 
ber deployment speed to M = .4 eliminated the controllability 
problems. 
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CHAPTER 5 


IMPLEMENTATION 

This chapter will describe the hardware and software mech- 
anization of the DAST I systems. The construction of the wing 
is described with details on the jettisonable mass ballast mechanism 
provided. The vehicle instrumentation is covered and the 
mechanizations of the FSS and RAV control systems described 


previously are detailed. 

5.1 ARW-1 Design and Construction . - The ARW-1 was de- 
signed by Teledyne Ryan Aeronautical of San Diego, California, 
under contract NAS1-13541 from NASA LaRC. 

The above contract specified the following design require- 
ments : 


The supercritical wing will be designed for cruise 
at a Mach number of 0.98 at 13,716 meters altitude 
and will have the following characteristics: 


Design lift coefficient 
Reference Area, Square Meters 
Thickness, Root, % chord 
Thickness, Tip, % chord 
Aspect Ratio 
Span , Meters 

Sweep at 25% chord, Degrees 
Sweep at 50% chord. Degrees 
Taper Ratio 

Design Load Limit at 1041 kgs 
(with a 1.25 factor of safety) 


.36 

2.787 

11.0 

7.0 

6.8 

4.343 

42.24 

40.0 

.36 

+2.5 

-1.5 


The outboard sections of ARW-1 are a 3/8 scale model of 
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the F-8 supercritical wing (reference 2) , but the wing-fuselage 
fairing, called the wing glove, is specifically adapted for 
the DAST. Access doors and spring loaded sway brace doors were 
designed for the glove section for equipment access and clear- 
ance for the forward sway braces on the wing pylon of the launch 
vehicle. 

The wing glove was designed of laminated fiberglass be- 
cause of the extreme contour changes. It is designed with 
aluminum formers and only minimal stresses are transferred to 
the glove/wing and glove/fuselage joint fasteners. 

The wing center box is machined from thick aluminum plate 
and was not optimized for weight. The spar material is steel 
because of the subscale size restrictions of the vehicle. The 
ribs are formed aluminum except the tip rib, which is machined 
from solid stock. "Tee" clips for the wing rib/spar joints are 
located where shallowness of the spar will not permit an inte- 
gral spar flange type of attachment. The leading and trailing 
edges are of laminated fiberglass construction to minimize 
torsional stiffness and for ease of manufacture. They are 
removeable for access to the wing instrumentation. The trail- 
ing edge is foam filled for skin panel support. The wing tip 
is a removable laminated foam filled fiberglass molding for 
instrumentation access. The laminated fiberglass skins be- 
tween spars are attached permanently with rivets. All detach- 
able portions of the wing and glove skins are attached by 
Voi-Sham screws . An illustration of the construction is shown 
in figure 5.1. 
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TOP SKIN 




FIGURE 5.1 


ILLUSTRATION OF WING CROSS-SECTION 



The wing was constructed at the NASA LaRC Model Shop 
in Hampton, Virginia. A picture of the completed wing is seen 
in figure 5.2. The leading and trailing edges are removed 
from the left wing to show the rib arrangement. 

The jettisonable mass ballast, which was added to lower 
the open loop flutter boundary when it was found that the 
fiberglass skin was stiffer than predicted, is shown in figure 
5.3. The lead shot in this cylinder is ejected when the re- 
taining cap is pushed back approximately 2.5 centimeters by a 
pyrotechnic activated bellows. The pyrotechnic charge is fired 
when an accelerometer signal exceeds a predetermined trigger 
level for a finite time interval. The electric circuit con- 
trolling the jettison mechanism integrates the accelerometer 
signal whenever it exceeds the trigger level of approximately 
20 g's. When this integrated value reaches a preset level, 
the pyrotechnic bellows are fired. The integration scheme pre- 
vents spikes in the accelerometer signal from firing the charge 
and insures that the wing is indeed fluttering before jettison- 
ing the shot. In tests on this system at NASA LaRC, the time 
required to fire the pyrotechnic charge once oscillations at 
the ARW-1 flutter frequency exceeded the trigger level was 0.4 
of a second or approximately six cycles. 

The GVT tests on ARW-1 and the DAST I vehicle arrangement, 
which was first mentioned in section 3.1.1, were conducted on 
the vehicle shown in figure 5.4. Only the forward wing glove/ 
fuselage fairings are missing from what represents the DAST I 
that will be tested at NASA DFRC. 
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FIGURE 5.2 COMPLETED ARW-1 
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FIGURE 5.3 THE JETTISONABLE TIP MASS BALLAST ATTACHMENT 
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FIGURE 5.4 DAST I GROUND VIBRATION TEST ARTICLE 
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5.2 Instrumentation . - The primary telemetry information 
for the real-time flutter monitoring during the flight tests 
will be from accelerometers in the wings and potentiometers 
measuring FSS control surface movements. The PCM main frame 
rate is 500 samples per second. This high sample rate was 
deemed necessary for the adequate time domain reconstruction 
of the higher frequency modes in the response. 

The rigid body sensor signals required for the RAV control 
system will be subcommutated at 100 and 50 samples per second. 

The data will be displayed on strip chart recorders in the 
DFRC control rooms . PCM tapes will also be generated for post 
flight analysis. 

Mach number will be computed from static and dynamic pressure 
measurements from the compensated probe attached to the nose 
cone. This probe is the same one that was used on the F-8 SCW 
flight tests. The nose boom is position error compensated to 
account for the DAST I f lowfield. Two Datametric pressure 
transducers were required to achieve the proper Mach number 
resolution in the transonic region. Resolution in Mach number 
will be at +.004 but absolute accuracy will only be achieved 
through dedicated airspeed calibration flights (which are not 
planned) . 

Ten accelerometers are involved in the FSS test. These 

are : 


Location 


Left Wing FSS Accelerometer 
Right Wing FSS Accelerometer 
Left Wing Front Spar Accel. 
Left Wing Rear Spar Accel. 


WBL-76 

WBL-76 

WBL-79 

WBL-79 
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Location 


Right Wing Front Spar Accel. 
Right Wing Rear Spar Accel. 

FSS Fuselage Normal Accel. 

FSS Fuselage Roll Angular Accel. 
Left Wing Tip Mass Accel. 

Right Wing Tip Mass Accel. 


WBL-79 
WBL-79 
FS 265 
FS 274 
WBL-79 
WBL-79 


The left and right wing aileron positions are measured by 


potentiometers attached to the actuator shafts. 

The rigid body sensors are illustrated in the on-board 
equipment layout in Section 5.4.1 except for the boom mounted 


vanes measuring angle of attack and sideslip angle. 

The pertinent instrumentation parameters and their accuracies 
are listed in Table 5.1. 

5.3 FSS Mechanization . - An illustration of the FSS layout 
on board the vehicle is shown in figure 5.5. The hydraulic 


servovalves were placed in the fuselage due to size constraints 


in the wing. The relatively long hydraulic lines to the act- 
uators near the wing tip were analyzed as part of the actuator 
math model and did not create a problem. The hydraulic pump, 
an electrically driven Sundstrand-Pesco Model 165-100, is rated 
at 3.558 x 10 -2 cubic meters of fluid per minute at 1.034 x 
10 7 Newtons/meter 2 . It was chosen because it met the FSS 
requirements and was the only "off the shelf" pump to do so 
that was compatable with the DAST I electrical power supply 
capabilities. This pump can supply 33.5% of the expected peak 
flow requirements. A 4.097 x 10 -4 cubic meter capacity accumu- 
lator will supply the remainder of the peak flow demand. 

The actuator location in the wing is shown in figure 5.6. 
This actuator is a single vane rotary design using an aluminum 
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TABLE 5 . 1 
DAST I 

TM PARAMETER 

SYSTEM I 
(PAGE 1 OF 3 


PARMID 


FSS 

Signals 

Name 

1. 

LWPOT 

LW Aileron Po.s 

2. 

RWPOT 

RW Aileron Pos 

3. 

FSSEXC 

FSS Excitation 

4. 

FSSLAC 

LW FSS acc. 

5. 

FSSRAC 

RW FSS acc. 

6. 

LWSERV 

LW servo command 

7. 

RWSERV 

RW servo command 

9. 

FSFUVAC 

FSS Fuselage A z 

9. 

FSFURAC 

FSS Roll ang. acc. 

10. 

DPAR 

FSS 'D' parameter 

11. 

HYDSUP 

FSS hyd supply press 

12. 

HYDLOD 

FSS hyd diff. load p: 

System II signals 


13. 

LFRACCL 

LW front spar acc. 

14. 

LRRACCL 

LW rear spar acc. 

15. 

RFRACCL 

RW front spar acc. 

16. 

RRRACCL 

RW rear spar acc. 

Tip 

mass signals 



17 . IjTMRAC 


LW tip mass acc. 


LIST 


) 


Range 

SPS 

Filter 

+12° 

500 

6P,70hz 

+12° 

500 

6P/ 7 Ohz 

T4° 

500 

6P,70hz 

+10g 

500 

3P , 70hz 

+10g 

500 

3P, 70hz 

+ 12° 

250 

3P,70hz 

+12° 

250 

3? , 7 Ohz 

+ 5g 

250 

3P,70hz 

+ 4 g 

250 

3P , 7 Ohz 

0-1000 units 

50 

IP , 40hz 

0-1500 psid 

50 

IP , 40hz 

+1500 psid 

50 

lp , 4 Ohz 


TBD 

500 

6P , 70hz 

TBD 

500 

6P, 70hz 

TBD 

500 

6P , 70hz 

TBD 

500 

6P,70hz 


250 3P,70hz 


+ 50g 
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TABLE 5.1 ( CONTINUED) 

FW tip mass acc. (PAGE 2 OF 3) +50g 250 


Avionics signals 


19. 

X001 

LM stab pos 

8,-13° 

50 

lP,40hz 

20. 

X002 

RH stab pos 

8,-13° 

50 

lP,40hz 

21. 

X003 

Rudder pos 

+10° 

50 

lP,40hz 

22. 

X004 

throttle pos 

0-270° 

50 

lP,40hz 

23. 

X005 

LH stab EAC 

8,-13° 

50 

IP, 40hz 

24. 

X006 

RH stab DAC 

8,-13° 

50 

lP,40hz 

25. 

X007 

Rudder DAC 

+10° 

50 

lP,40hz 

26. 

X008 

Throttle DAC 

0-270° 

50 

lP,40hz 

27. 

M001 

Pitch angle 

+80° 

50 

lP,40hz 

28. 

MO 02 

Roll angle 

+] 80° 

50 

lP,40hz 

29. 

M003 

yaw angle 

+] 80° 

50 

lP,40hz 

30. 

M004 

RPM 

0-110% 

50 

lP,40hz 

31. 

M005 

Fuel flow rate 

0.5/8 gpm 

50 

lP,40hz 

32. 

M006 

Fuel remaining 

0-305 lbs. 

50 

lP,40hz 

33. 

T001 

EGT 

700-850°C 

50 

lP,40hz 

34. 

M007 

AGC 

0/-110db 

50 

lP,40hz 

35. 

MO 08 

Alpha 

+30,-12° 

100 

lP,40hz 

36. 

M010 

Beta 

+20° 

100 

lP,40hz 

37. 

M011 

+12 VDC monitor 

0-25VDC 

50 

lP,40hz 


Instrumentation Signals 


38. 

A0 01 

Pitch rate 

+60 °/s 

100 

lP,40hz 

39. 

A002 

roll rate 

+100 °/s 

100 

lP,40hz 

40. 

A003 

yaw rate 

+40 o/ s 

100 

lP,40hz 

41. 

M012 

+28 VDC monitor 

0-30 VDC 

50 

lP,40hz 

42. 

M013 

26 VAC monitor 

0-28 VAC 

50 

lP,40hz 

43. 

A004 

Az,cg 

+6,-3g 

100 

lP,40hz 

44. 

A005 

Az,cg fine 

+3,-0. 5g 

100 

lP,40hz 

45. 

A006 

Ay, eg 

+1 g 

100 

lP,40hz 

46. 

A007 

Ax, eg 

+i g 

100 

lP,40hz 

47. 

M009 

Alpha, fine 

+6,-2° 

100 

lP,40hz 

48. 

D001 

Airspeed diff. press. 

0-1500 psfa 

50 

lF,40hz 
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49. STADIF 


50 


lP,40hz 


TABLE 5.1 (CONCLUDED) (PAGE 3 OF 3) 


Ps-Pref , AP 


+2.5 psid 


Digital Signals - Avionics 


1 . 

M014 

Recovery 

2. 

M015 

AP 

3. 

PHOCCM 

Camera on-off 

4. 

FSSON 

FSS on-off 

5. 

PUICCM 

FSS pulse excitation 

6. 

M0 16 

Emergency Recovery 

7. 

1MRC0M 

Tip Mass eject. 

8. 

M0 17 

smoke 

9. 

AMPOOM 

FSS excitation amp. 

10. 

SWEOOM 

FSS sweep excitation 


DW-1 

DW-1 

DW-1 

DW-1 

DW-1 

DW-1 

DW-1 

DW-1 

DW-1 

DW-1 


50 

50 

50 

50 

50 

50 

50 

50 

50 

50 


Digital Signals - Avionics 


1 . 

FSGAIN 

FSS gain 

2. 

SYMCCM 

FSS excitation 

3. 

T .KkKl i 

IK Pyro verify 

4. 

RIREL 

EN Pyro verify 

5. 

HRESLV 

Hyd pump Res. status 

6. 

M018 

90° roll switch 

7. 

PMPON 

Hyd pump on-off 


DW-2 hi-low 50 
DW-2 sym-anti 50 
DW-2 50 
DW-2 50 
DW-2 50 
DW-2 50 
DW-2 50 


Digital Signals - Inst ♦ 

1. P002C altitude, course (6kits) 

2. P002F altitude, fine (6 kits) 


EW 3 
DW 4 


50 

50 


PCM Frame Status - Channels 


Available 


Used 


500 SPS 
250 SPS 
100 SPS 
50 SPS 


9 

6 

10 

30 


9 

6 

10 

28 (24 analog, 4 digital) 


123 



0.907 KG 
0.15 M 


FIGURE 5.5 


FSS ON-BOARD LAYOUT 
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ACTUATOR 


I 



FIGURE 5.6 


CONTROL SURFACE AND ACTUATOR ILLUSTRATION 


The control surface 


body with a stainless steel shaft and vane, 
is directly coupled to the shaft. This system was sized so that 
the stall torque equals the maximum expected hinge moment. The 
actuator has a 5.03 meter/sec peak velocity at 1460 deg/sec. 

The servoactuator block diagram, shown in figure 5.7, shows the 
closed loop transfer function which was used in the analysis 
in Chapter 6. 

The functional interface requirements of the FSS with the 
rest of the DAST I systems is illustrated in figure 5.8. Dual 
potentiometers are utilized on the actuator shaft so that a 
spare would exist. A possible major disassembly of the system 
is prevented by incorporating the spare into the actuator at 
the initial buildup. The two pressure tranducers are utilized 
with the servovalves so that absolute and differential pressure 
can be monitored. 

5.4 Flight Control Mechanization . - The flight test philo- 
sophy for the FSS experiment is based on the vehicle operating 
under RAV control. If the telemetry link required by the RAV 
control system is broken, automatic reversion to the on-board 
autopilot occurs. This interface is critical to safe operation 
of the DAST RPRV. 

5 ^ 4.1 PAST I On-Board Systems . - One of the critical systems 
in the operation of the DAST as an RPRV is the wrap-around nose 
cone antenna. The basic Firebee drone had a receiving and trans- 
mitting antenna on the lower surface of each wing tip. The single 
nose cone antenna for the DAST I serves as both the point of 
transmission and receiving. The nose cone antenna was fabricated 


125 


126 


SERVOVALVE ACTUATOR, LINES AND CONTROL SURFACE 



CLOSED LOOP NO-LOAD TRANSFER FUNCTION 

2.721 x IQ 21 (S +10) DEG/DEG 

<* c (S+ 10.14) (S + 375 . 4) (S 2 + 758. OS + 7.158 x 10 5 )(S 2 + 852. 9S+ 1 .591 x 10 6 )(S 2 + 1164S+ 6.277 x 10 6 ) 


FIGURE 5.7 


FSS SERVOACTUATOR BLOCK DIAGRAM 
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FIGURE 5.0 


ARW-1 FSS INTERFACES 








at NASA LaRC and tests were performed to determine its signal 
strength pattern. This pattern is illustrated in figure 5.9. 
The approximate radius of sufficient signal strength of the 
pattern is on the order of 96.5 kilometers. The furthest 
planned distance from the ground transmitting and receiving 
station is 72.4 kilometers. The "holes" in the pattern which 
are evident forward and aft of the vehicle, are potential 
sources of the TM "dropouts" as the vehicle turns toward and 
away from the ground station. This likelihood increases as 
the vehicle's altitude decreases. 

Space available in the fuselage for the vehicle systems 
has been a constraining factor. The on-board equipment layout 
is shown in figure 5.10. The brackets, wiring, and plumping, 
which are not represented here, make access to the on-board 
equipment difficult so that system removal for debugging and 
repair is a lengthy process. 

The on-board control system mechanization and RAV/AP in- 
terfacing is illustrated in figure 5.11. The 10-bit digital 
proportional words described in section 5.5.3 are transmitted 
by the RPRV ground facility or the F-104 secondary command 
station. If the autopilot has to be activated, either auto- 
matically by a loss of telemetry or by manual switching, then 
the commands are processed by the autopilot. If the system 
is in the RAV mode, then the autopilot is disabled. When a 
switch to the autopilot occurs, the RAV commands are ramped 
out over a three second interval to minimize the transients 

Recall that after the RAV/AP switch, the 


in the pitch axis. 
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FIGURE 5.9 DAST I NOSE CONE ANTENNA PATTERN 
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FIGURE 5.10 


ON-BOARD EQUIPMENT LAYOUT 
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FIGURE 5.11 DAST 0N-30ARD CONTROL SYSTEM INTERFACE 


autopilot will drive the system to a zero glideslope condition 
without F-104 or ground station inputs. The discrete words, 
which will be described in section 5.5.3 and are also uplinked 
to the vehicle, control the functions identified in the sche- 
matic shown in figure 5.12. The uplink commands consist of 
four 16-bit digital words. The first 10 bits of each word con- 
tain the proportional commands. The remaining six bits on each 
word (24 total) are packed with discrete commands. The on- 
board data interface box separates the proportional and discrete 
commands and routes them to the appropriate locations. The 
discrete commands for smoke, recovery, and emergency recovery 
can be changed by internal timers which are activated if an 
automatic transfer to the AP mode occurs. If commands to the 
DAST are not received from the F-104 chase aircraft within 
the 20 second limit of these timers, the automatic recovery 
begins . 

5.4.2 F-104 Controller . - The F-104 secondary command 

station will only take control of the DAST I if the telemetry 
link with the RPRV facility on the ground is lost. The fail- 
ure in the telemetry link must be in the ground station or be 
an on-board downlink failure, otherwise the F-104 cannot con- 
trol the DAST either. This is obvious when one considers that 
the downlink telemetry is not required for F-104 control, 
whereas reception of the uplink signal by the vehicle is re- 
quired for control by both the RAV and autopilot systems. If 
the TM failure is due to a reception problem on-board the RPRV, 
automatic recovery will result. 


132 



FIGURE 5.12 


ON-BOARD TM PROCESSING AND DISCRETE COMMANDS 


The RPRV pilot in the rear seat of the F-104 can control 
the DAST engine rpm and make climb, dive, and bank angle com- 
mand inputs to the vehicle. Figure 5.13 shows the RPRV control- 
panel in the F-104. There is an automatic throttle reduction 
to the 80% rpm setting whenever the AP mode is activated. If 
the RPRV pilot determines that more thrust is necessary , the 
dial on the panel is lifted and rotated. The F-104 chase air- 
craft RPRV pilot can only judge the throttle requirement re- 
lative to the F-104 airspeed since he will have no indication 
of the magnitude of engine rpm changes he makes . 

Turns and altitude changes are made by the trim switch on 
the panel. The integration rates and authority limits on this 
stick are such that no high g maneuvers can be made and the 
vertical velocity is limited to a maximum of + 30 meters per 
second. The bank angle command is limited to + 40 degrees. 

A straight and level switch on the panel, when utilized, will 
zqto the pitch and roll commands and let the autopilot main- 
tain altitude and heading. Extensive piloted simulation was 
performed to verify adequate and safe values for the trim rates 
and limits throughout the flight envelope. 

5.5 RPRV Ground Facility . - The NASA DFRC RPRV facility 
was developed to enable high risk flight testing. Due to the 
nature of some flight experiments, subscale remotely piloted 
aircraft were chosen as the safest and most efficient method 
of performing some of these tests. This section will describe 

this facility and how it operates. 

5 . 5.1 Ground Cockpit . - The ground cockpit used for the 
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FIGURE 5.13 F-104 RPRV CONTROL PANEL 
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the basic Firebee RPRV is shown in figure 5.14. This is the 
cockpit which will be used with the DAST I with some changes 
in instrument displays. It is a conventional center stick and 
rudder pedal configuration. The forces and displacements 
chosen for the DAST I pilot commands are shown in figure 5.15. 
Electric motors controlled by the cockpit stick computer, gen- 
erate the artificial feel system forces shown in figure 5.15. 
These gearing schedules are the same as those that were utilized 
on the Firebee RPRV. 

The cockpit instrument displays and radar tracking plot 
board with an area map overlay combine to form the information 
the pilot will use to fly the DAST. The information displayed 
to the pilot from the cockpit instruments are: 

1. angle of attack 

2. sideslip angle 

3. Mach number 

4. airspeed 

5. 3-axis gyro displaying pitch attitude, roll attitude 
and heading 

6. altitude 

7. vertical velocity 

8. normal acceleration 

9. surface positions on strip gages 

10 . surface commands 

11. engine rpm in percent 

12 . commanded engine rpm 

13. fuel flow rate 

14. fuel quantity in weight 

15. engine exhaust temperature 
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The fuel quantity is determined from integrating the fuel 
flow rate on board the vehicle. Mach number, airspeed, and 
vertical velocity are computed in the RAV computer from static 
and dynamic pressure measurements in the TM downlink. 

The mode control box seen to the left of the cockpit in 
figure 5.14, is the means by which the pilot can change gains 
and RAV modes during flight. Figure 5.16 is a more detailed 
illustration of the mode control box. Different modes are 
activated by depressing the appropriate button which lights 
up to verify that the RAV computer received the mode change 
and is executing the proper control law software. Gain sel- 
ection switches are seen at the bottom of each row. They also 
have lights which verify that the gain selected by the pilot 
is indeed the one being used by the RAV computer. 

All the lights on the mode control box are controlled by 
the control law software. Discretes from the buttons and gain 
switches are used to activate certain segments of software 
code. The box also allows selection of a computer bypass mode 
and informs the pilot if any downlink variables fail a window 
check. 

The pilot has numerous switches in the cockpit which 
change the discrete words in the uplink to the vehicle. These 
are : 

1. RAV/AP Transfer 

2 . Recovery command 

3. Emergency Recovery 

4. Ground Transmitter On/Off 

5. Smoke Generation 
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FIGURE 5.16 RAV MODE CONTROL BOX 
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6. Control Surface Pulses and Doublets 

The ground transmitter on/off switch is used to turn the 
ground transmitter off when the F-104 is commanding the vehicle. 
The smoke generation function switch causes oil to be injected 
into the engine exhaust for enhancement of visible tracking. 

The pulse panel would be used to generate rudder, rolling tail, 
and elevator preprogrammed inputs of selectable magnitude for 
stability and control maneuvers. 

5.5.2 Ground Computer . - A photograph of the RAV ground 
computing equipment is shown in figure 5.17. This facility 
is the "heart" of the RAV concept. The control laws and other 
functions to be described in section 5.5.4 are executed in 
real-time. There are three interface systems linking the 
ground computer into the closed loop RAV operation. 

Telemetry Downlink Interface . - Downlink telemetry data 
are received at 500 samples per second. The variables re- 
quired for the computer are subcommutated and placed in a 
hardware storage buffer by a telemetry downlink computer. This 
information is transferred from the hardware storage buffer 
to the RAV computer on an as needed basis by a call statement 
in the control law software. The values of the variables are 
represented by a 0 to 511 decimal count format. 

During real-time operation, the downlink interrupt sync 
pulse, which occurs once every PCM frame, drives the downlink 
computer service routines to continuously update the downlink 
storage buffer so that the most recent data will always be 
available to the control laws computation. This minimizes the 
time delay in the closed loop cycle. 
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Telemetry Uplink Interface . - The computer encoder uplink 
interrupts the computer every 3.75 milliseconds to request 
one of four 16-bit data command words or a 16-bit synchroniza- 
tion word. If the computer program and uplink encoder get 
out of synchronization, the computer sends the variable that 
the encoder requested. The uplink mainframe rate is 53.33 
samples per second. 

Ground Cockpit Interface . - The pilot's proportional com- 
mand signals, mode control box, and cockpit switches status 
are sampled once each cycle at the 53.33 samples per second 
rate. The RAV software main FORTRAN program determines if 
any mode or gain changes are being commanded by the pilot. 

The cockpit data is transferred to the RAV computer during 
the synchronization word computation mentioned above. 

The computer used in the RPRV system is a general purpose 
rack-mounted mini-computer with a 32K memory of 16-bit words 
and with a 330-nanosecond cycle time. Peripheral equipment 
includes a CRT, card reader, line printer, magnetic tape unit, 
disc unit, and peripheral floating point processor. The soft- 
ware is composed of an assembler, a FORTRAN compiler, and a 
mathematical subroutine support library. The open-loop and 
closed- loop control law computations are implemented in 
floating point FORTRAN. This allows the FORTRAN compiler to 
be used in debugging and checking out programs , and the float- 
ing point feature eliminates the need for variable scaling. 

The FORTRAN feature is ideally suited for research programs 
because of the ease with which programs can be written and 
changed. 


144 


5.5.3 RAV Software Structure. - Because of the experience 


gained at NASA DFRC in the operation of RAV and other flight 
programs requiring flight software (F-8-DFBW aircraft) , rather 
stringent and extensive software control requirements have 
been established. For safety of flight considerations, ex- 
tensive verification testing and documentation procedures are 
required and reviewed by a designated program Software Control 
Board. 

This section is intended to provide a description of the 
required software structure which makes the RAV concept work- 
able while sections 5.5.4 and 5.5.5 will describe some of the 
safety features and possible failure areas which can occur 
from RPRV operations. 

The RAV program executive software structure is illustrated 
in figure 5.18. The background displays and program initial- 
izations are executed only in the non-real time mode of the 
program. 

The CRT background displays will be used for: 

• Internal program monitoring 

• Software and hardware debugging 

• Program status and timing information 

These diagnoistic and informational displays will be one 
of the primary tools used in the software verification tests. 
Line printer copies of these displays are available for doc- 
umentation. All the information required for proper RPRV fac- 
ility checkout and software functional checks will be avail- 
able on these CRT displays. 
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FIGURE 5.18 RAV PROGRAM EXECUTIVE STRUCTURE (SIMPLIFIED) 


146 


The program initiation will be performed prior to activa- 
tion of real-time operations. Its purpose is to ensure that 
no program variable or hardware element is left undefined 
prior to entering the real-time software code. Software 
variables which require initialization include limits, gains, 
constants, calibrations, and the fault thresholds, tolerances, 
and documentation of fault occurrences which will be discussed 
in section 5.5.4. The hardware elements which require initia- 
lization are digital- to-analog converters, analog-to-digital 
converters, and the uplink controller. 

The real-time functions of the RAV software are the real 
time loop structure, the system input/output interfaces, and 
the control laws. The Babcock uplink encoder external interrupt 
drives the real-time loop structure of the DAST computer program. 
This interrupt is on an 18.75 millisecond frame rate. During 
the uplink interrupt loop, the control law algorithms, some 
elements of the system input and output interfaces, and the fault 
detection algorithms are executed. The downlink operation 
requires no program external interrupts, although a fault de- 
tection algorithm will utilize the PCM sync word arriving 
every 2 milliseconds. 

The downlink interface described in section 5.5.2 between 
the RAV computer and the downlink servicing computer transfers 
the variables required in the computations in decimal counts, 
then converts them to engineering units using linear slope/ 
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intercept formulas before they are used in the uplink interrupt 
cycle. 

The downlink proportional inputs required in the computa- 
tions are: 

1. Pitch rate 

2. Yaw rate 

3. Normal acceleration 

4. Free stream static pressure 

5. Free stream differential pressure 

Several downlink discrete inputs are required by the ground 
computer to provide status information to the pilot. These 
are: 


1. 

Carrier loss discrete 


2. 

FSS hydraulic pump reserve status 

3. 

Left wing pyrotechnic verify 

(fired) 

4. 

Right wing pyrotechnic verify 

(fired) 

5. 

FSS gain verify 


6 . 

Recovery commanded 


7. 

FSS on/off verify 


8. 

Emergency recovery 



The four uplink data words generated are shown in figure 
5.19. The first ten bits of each uplink word are under com- 
puter control and the last six bits are hardwired to the up- 
link encoder. It can be seen from figure 5.19 that the dis- 
crete words are complimented on every other digital word. It 
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BITS 0 THRU 9 

BIT 10 

BIT 11 

BIT 12 

BIT 13 

BIT 14 

BIT 15 


RECOVERY 

AP CMD 


ANTI- 

CAMERA ON 

FSS ON 

PULSE ON 

WORD 1 LEFT HAND STAB 



SYMMETRIC 

0 (OFF) 

1 (OFF) 

0 (OFF) 

1 (OFF) 

0 (OFF) 

1 

(OFF) 


3MEH 

TIP MASS 



SMOKE 

AMPLITUDE 

FREQ ON 

WORD 2 RIGHT HAND STAB 

RECOVERY 

EJECTION 



0 (OFF) 

HI 

0 (OFF) 

1 (OFF) 

0 (OFF) 

l_L 

(OFF) 

t (OFF) 


RECOVERY 

AP CMD | 


ANTI- 

CAMERA ON 

FSS ON 

PULSE ON 

WORD 3 YAW 



SYMMETRIC 


0 (OFF) 

1 (OFF) 

0 (OFF) 

1 (OFF) 

0 

(OFF) 

1 (OFF) 


EMSR 

TIP MASS 



SMOKE 

AMPLITUDE 

FREQ ON 

WORD 4 THROTTLE 

RECOVERY 

EJECTION 



1 (OFF) 

HI 

1 (OFF) 

0 (OFF) 

1 (OFF) 

0 

(OFF) 

0 (OFF) 


FIGURE 5.19 UPLINK WORD STRUCTURE 




has been found that this scheme minimizes the changes of a 
spike in the uplink generating an unwanted command. 

The control law computation structure within the uplink 
interrupt cycle is illustrated in figure 5,20. Since the 
RAV control laws were not complicated, the computation time 
for the control laws is only requiring approximately 42 
per cent of the 18.75 milliseconds available. The com- 
putations for driving the Mach meter, airspeed indicator, 
vertical velocity indicator and fault detection algorithms 
therefore, do not tax the system. These calculations are 
divided up so that their software call statements are evenly 
distributed with respect to the control surface computations 
of the control laws. 
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LOOP 



FIGURE 5,20 CONTROL SYSTEM EXECUTE STRUCTURE (SIMPLIFIED) 
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5.5.4 Failure Detection Methods . - Although the RPRV con- 
cept has proven to he an excellent method of performing high 
risk flight experiments, the concept itself introduces addi- 
tional mechanisms from which a failure in flight can occur. 
Fortunately, the computational capability which is available 
in the ground facility can reduce the chances of endangering 
the vehicle from some failures which can occur with the RPRV 
technique . 

The mechanisms which have been employed in past and pre- 
sent RAV programs at NASA DFRC and which will be used in the 
DAST program are summarized in Table 5.2. The time-out 
counter concept is illustrated in figure 5.21. The purpose 
of the time-out counters are to prevent asynchronous opera- 
tion of the uplink and downlink during closed-loop operation. 
Asynchronous operation would allow unacceptable time delays 
in the control laws and other possible dangerous operations. 

The rate checks on the uplink commands prevents "hard 
over" commands, that could result from something like a 
short circuit in the cockpit stick sensors , from being sent 
to the vehicle. Downlink rate checks protect against possible 
sensor "hard-over" failures from being used in the feedback 
computations. Automatic mode changes and AP transfer are 
ways of preventing dangerous and unintentional RAV commands 
resulting from the types of failures described above from 
being sent to the vehicle. 

5.5.5 Failure Modes Assessment - Part of the policy of 
qualifying software for flight is to perform a failure modes 
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TABLE 5.2 

RAV Fault Detection Mechanisms 


Mechanisms 

Detects 

Downlink Software Time-Out Counter 

Hang Up in Downlink Process 

Uplink Software Time-Out Counter 

Hang Up in Uplink Process 

Uplink Command Rate Check 

Unreasonable Changes in RAV Commands 

Uplink Command Authority Limit 

Out of Range Commands 

Downlink Variable Rate Check 

Unreasonable Changes in Downlink 
Parameters 

Stick-Trim Opposition Logic 

Trim Discrete Failure 

Computer Heartbeat Monitor 

Ground Computer Hardware Malfunction 

Pilot 

Unexpected, Uncommanded Vehicle Motion 









and effects analysis (FMEA) on the system. Prior to be- 
ginning a FMEA, the software must have successfully completed 
static, dynamic, and closed- loop real-time simulation testing. 
Static checks are accomplished by holding all inputs constant 
and checking the output of each control system element for 
proper values. In dynamic verification testing, all dynamic 
control system elements, such as filters, rate checks, and 
faders, are checked individually input to output. Closed loop 
operation is verified in the real-time simulation. When the 
software has been verified to this level and documented, an 
FMEA is the next step. The proper operation of the failure 
detection schemes are verified in the types of test shown 
in Table 5.3. 

An on— board receiver failure, autopilot sensor failure, or 
horizontal stabilizer servo failure all place the vehicle in 
potentially disasterous circumstances. The automatic re- 
covery resulting from an on-board receiver failure is consid- 
ered potentially disasterous because it will probably mean the 
recovery parachutes will deploy outside the recovery area. 

This would mean that the MARS helicopter could not reach the 
DAST before it impacted the ground. This is also true if a 
horizontal tail servo fails and the vehicle survives any un- 
controllable maneuvers before the recovery chutes deploy. With 
an autopilot sensor failure the vehicle would probably en- 
counter uncontrollable maneuvers when the vehicle recovery is 
commanded, thus jepardizing the success of the recovery. 
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TABLE 5.3 

Failure Modes and Effects Analysis 


Simulation 


RPRV Facility/DAST 


Downlink Faults 


• opens 

• hardovers 

• ramps 

• dropout 

• frozen value 

Uplink Faults 


• 

open 

• 

hardovers 

• 

ramps 

# 

frozen value 

Software Faults 

• 

hardover 

• 

ramps 

• 

frozen value 

• 

runaway integrator 

| 

Ground Computer 

# 

power loss 

• 

halt 


Downlink Faults 


• opens 

• hardovers 

• ramps 

• frozen value 

• TM station 

• power loss 

• TM station 
sync loss 

Uplink Faults 


• opens 

• hardovers 

• ramps 

• frozen value 

• signal loss 

• parity fail 

• power loss 

• Babcock encoder 
failure 

• TM transmitter 
failure 
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CHAPTER 6 


ASSESSMENT OF FLUTTER TESTING METHODS 

This chapter examines some of the unique aspects of the 
FSS flight testing requirements. A testing philosophy is 
discussed which addresses the new aspects of flutter testing 
this research experiment creates. An evaluation of candidate 
testing techniques is accomplished by using DAST ARW-1 sim- 
ulated flutter data. The development of the DAST ARW-1 
flutter simulation is also described. 

6.1 Flight Testing Philosophy and Limitations . - The 
FSS flutter testing experiment is a new regime of flutter 
testing. The flutter characteristics are explosive and of 
high frequency. Demonstrating an increase in the flutter 
speed of this magnitude has never been attempted before in 
flight testing. For DAST I, the testing philosophy is to 
rely on the tip mass ejection system to prevent the wing 
from fluttering if problems arise in the supercritical test- 
ing region. A test engineer on the ground will monitor the 
accelerometer signals but since he will not be able to identi 
fy a flutter problem on this wing in sufficient time to man- 
ually jettison the tip ballast, the automatic system will be 
relied upon to catch any problems. Only closed loop testing 
of DAST I will occur near the flutter boundary. The risk, 
because of the explosive flutter, makes precise determination 
of the unaugmented boundary very dangerous. It will be de- 
termined "as well as possible" so that some gage of the FSS 
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performance can be made. A preliminary first flight test 
point grid is shown in figure 6.1. The lowest altitude test 
points will be completed first so the vehicle can decelerate 
while climbing to the next test altitude. 

A possible second flight test grid is seen in figure 
6.2. The flight plan calls for the FSS to be on at all times. 
The unaugmented flutter boundary will be exceeded on the sec- 
ond flight. Consideration is being given to demonstrating 
the 20 per cent increase in the flutter boundary on the third 
flight. The test program is scheduled for six flights at a 
one flight per month rate. 

The difficulty of conducting the DAST I FSS experiment 
is increased slightly by limitations introduced by the ve- 
hicle's performance characteristics. These are: 

1. Limited flight time 

2. Limited antenna range 

3 . Lack of precise speed control 

4. Limited FSS hydraulic capacity 

5. The relatively low MARS reliability 

In examining the first three items listed above, their 
effects on the flight tests can be minimized by adopting 
flutter testing techniques which require short record lengths 
of data. 

The limited flight time, caused by a small fuel capacity 
and a large fuel consumption rate, is still of concern but 
the amount of data per flight is greater with these new test- 
ing methods. With the limited antenna range and high vehicle 
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FIGURE 6,1 FIRST FLIGHT TEST POINT GRID 
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FIGURE 6.2 SECOND SUPERCRITICAL FLIGHT TEST POINT GRID 
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speeds, the DAST can fly only three or four minutes before a 
turn must be made. Since the data record length requirements 
are short, the constant flight condition requirements during 
data collection are easily achieved. The precise speed con- 
trol problem is resolved somewhat because the testing methods 
require data records on the order of only ten seconds . Very 
little change in Mach number will occur in this short time 
interval even if the vehicle has not been totally speed stab- 
ilized by the pilot. Of course, this is only true if the 
RPRV pilot is trying to establish nonaccelerating conditions. 
The point is that the flutter testing methods don't require 
the vehicle to be totally speed stabilized during the tests. 
Trying to establish a totally constant Mach number condition 
with DAST I is difficult because of the poor engine response 
and the rate of change in flight parameters due to the rapid 
fuel burnoff. 

The concern of limited FSS hydraulic capacity is created 
by the fact that 80 per cent of the FSS hydraulic capacity is 
required just to maintain the control surface in the faired 
position. Aft loading from the trailing edge camber of the 
SCW is the cause. If turbulence is encountered and excitation 
inputs are also added, the hydraulic system could be over- 
powered and surface saturation could result. Flying in rea- 
sonably calm flight conditions reduces the risk of this hap- 
pening. Another concern is that the load on the control sur- 
face in the faired position requires nonlinear actuator loop 
gain. This is illustrated in figure 6.3. This type of non- 
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FSS NONLINEAR ACTUATOR LOOP GAIN 



linearity introduces the possibility of a limit cycle in the 
FSS system. 

Air Force experience with MARS recoveries of Firebee II s 
shows an approximately 80 per cent success rate. This justifies 
the rather bold approach to the flight testing. By having 
a success oriented approach, the amount of data on the FSS 
performance characteristics is maximized. If, on the other 
hand, a more conservative and cautious philosophy were adopted, 
the chances diminish of completing the tests successfully 
before a MARS failure. By using some of the relatively new 
flutter testing techniques described in the next section, a 
large amount of information on the FSS system can be attained 
quickly. Thus, the risks involved in rapid expansion of the 

DAST flutter envelop are reduced. 

6.2 Near Real-Time Testing Methods . The requirement 
for flutter testing methods which can make estimates of 
structural frequency and damping from short record lengths 
of data has been established. An evaluation of some of the 
relatively new methods was done in reference 11 . Their per- 
formance was measured by testing the accuracy of frequency 
and damping extimates in a two-dimensional flutter model. 

This section concerns itself with describing those techniques 
which are suited to the requirements of the DAST experiment. 

Most flutter test frequency and damping estimation methods 
operate in the frequency domain. The frequency response trans- 
fer function, Ri, of an accelerometer signal is determined by 
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taking the Fast Fourier Transform (FFT) of the time domain 
acceleration response. 

Those methods considered as candidates for the DAST near 
real-time requirement are: 

1. PSD's 

2 . Co-Quad 

3. Amplitude-Phase 

4. Modulus of the Transfer Function 

5. Impulse Response Function 

The Power Spectral Density (PSD) technique is illustrated 
in figure 6.4. Very short data record lengths can be used 
with this method. The Spectrum of H 2 , which assumes white 
noise input , shows the modal frequencies by peaks in the 
power spectrum. The damping estimate is made by determining 
the width of the frequency peak at its half power point and 
dividing by twice the peak frequency. The white noise input 
is usually approximated by pulses and doublets of the control 
surfaces which create broad band width excitation. 

The Co-Quad method generally uses a frequency sweep of 
the control surfaces to excite the different modes. As the 
sweep frequency coincides with a modal frequency, energy is 
input into that mode and its response increases. The transfer 
function, H 2 is determined for this input and its real (coin- 
cidence) and imaginary (quadrature) components are plotted as 
in figure 6.5. When a structural mode is excited by the in- 
put, the amplitude of the transfer function increases. Then 
the damping estimate is determined from the changes in the 
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imaginary part at the half power frequencies. 

The PSD and Co-Quad methods assume that the form of the 
damping is structural. This means that the damping is pro- 
portional to the square of the amplitude of displacement of 
the structure. If viscous damping is present, the damping 
estimates produced from the PSD plot and Co-Quad plot will 
not be truly representative of damping within the system. 

The Amplitude-Phase method is not dependent upon the 
structural damping assumption. Figure 6.6 illustrates this 
method. The frequency response transfer function is used in 
its polar coordinate form. The amplitude ratio indicates 
the modal frequency and the damping is determined from the 
slope of the phase angle at that frequency. 

The modulus of the frequency response function, e.i., 
the square of the transfer function, will generate PSD's 
from system inputs such as surface frequency sweeps. Fre 
quency and damping are then determined in the same manner as 
the PSD method. Its advantage over the PSD method is its 
noise reduction capability which is discussed later in this 

section. 

The impulse response function is the inverse Fourier 
transform of the frequency response function: 


h (t) 


1 _ 

2it 




( 6 . 1 ) 
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Frequency and damping measurements from the impulse 
response are illustrated in figure 6.7. Frequency is measured 
by the period of oscillation and damping is determined from 
the log-decrement of the oscillation. 

Aside from the PSD method, all the methods deal with the 
noise problem by the following mechanism. The response due to 
an applied excitation and due to a noise environment may be 
written as: ~ 

Y = Y + Y = I (f + f_) *h (t-x)dT (6.2) 

f n J n 

where f is the unknown input noise and Y r is the associated 

♦ • 

noise in the response. The Fourier transform of Y is: 

Ft* = Ft: + Fit = (Ff + F-f ) H 2 (6.3) 

Y Y f n n 

If F is multiplied by the complex conjugate F f , then the 

f 

following spectrum equation results: 

^F” + ^Fit = (<f> + <t> ff ) h 2 

Y f x n f 1 n 

Since f and f are uncorrelated, the cross spectra, 

6 ff n , should vanish. This yeilds: 

't’p.. 

H 2 - Y f 

tp ^ 

which should appear as a noise-free result. 

For the analysis performed in Section 6.3, the frequency 


(6.4) 


<j) and 
Y n 


(6.5) 
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ACCELERATION IN G'S 



FIGURE 6.7 ILLUSTRATION OF THE IMPULSE RESPONE FUNCTION 


response function, H 2 , was generated in this manner. The 
impulse response function was determined from the inverse 
Fourier transfer of H 2 . The PSD's determined in the next 
section, are from simulated accelerometer responses only. 

For the DAST I flight tests, a modal analysis 
computing facility will be used to calculate the frequency 
response function and will be able to produce plots of any one 
of the methods described previously in about 10 seconds. 

Strip chart readouts of accelerometer signals will be used 
to watch for "ballooning” in the response. The communication 
link between the modal analysis facility, the control room, 
and the pilot in the RPRV facility is critical. The indivi- 
duals in all three locations must be well coordinated m their 

activities. 

Flutter margins, which all three groups must be conscious 
of, will be displayed on a plot board showing o> and S versus 
attitude and Mach number. Frequency coalesence will be evi- 
dent in this type of display and projections of the damping, 

S, to the zero crossing can indicate the altitude and Mach 
number that the flutter boundary might first be encountered. 

6.3 simulated Responses . - In an attempt to evaluate 
the flutter testing methods to determine which are best 
suited for the DAST experiment, a simulation of the DAST 
ARW-1 structural and rigid body motions was developed. The 
math model of the DAST ARW-1 developed for the FFS analysis 
was transformed into a state space representation and im- 
plemented in a modified version of the CONTROL computer 
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program. The transformation equations required to make this 
conversion are presented in Appendix C. 

Many of the data characteristics which the DAST will 
possess were duplicated in the model. These are: 

• 500 samples per second instrumentation 

• Identical control inputs 

• Record lengths not exceeding 4096 points 

• Noise (Dryden Turbulence) 

The unaugmented equations of motion were of 54th order. 

They include: 

• 1 rigid body mode 

• 6 wing and fuselage structural modes 

• 36 artifical states representing the unsteady 
aerodynamics 

• Gust effects 

• Control surface mode 

The augmented equations which include a simplified 
FSS control law, supplied by BWC, are of 60th order. The 
simplified FSS control law has characteristics very close to 
the complete FSS system up to 400 radians in frequency. 

The rigid body pitch mode was included to see if 
there was any appreciable coupling of the rigid body modes 
and the flexible modes. 

The data generated for the doublet and sweep inputs 
of the control surface are of the same record length as 
planned in the flight tests. The modal analysis system 
can only store and use a maximum of 4096 data points per 
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calculation. The record length of data points for the doublets 
is 904 to allow complete decay. The maximum 4096 points were 
used for the frequency sweep inputs. This required a time 
history of 5000 data points or seconds. 

The control surface inputs planned for the flight tests 
and used in the flutter simulation are a 25-hertz square wave 
doublet and a logarithmic frequency sweep from 10 to 40 hertz 
in 6 seconds. The sweep is ramped in and out over a .3 sec- 
ond time interval. These inputs are shown in figure 6.8. 

This is a time history plot of control surface commands and 
movements during a simulation run. PSD's of the surface 
movements in figure 6 . 8 show the frequency content of the in- 
puts. The 25-hertz doublet creates a broad bandwith input 
about the doublet frequency. The PSD of the control surface 
sweep shows the effect of the ramping in and out of the sweep. 

The power spectrum indicates that ramping in and out reduces 
the bandwith of the original input. If this were not done, however 
and the commands were inserted abruptly and removed abruptly, 
the spectrum characteristics at the beginning and end of the 
sweep would appear as large bandwidth peaks since the abrupt 
changes are similar to pulses. 

With the simulated responses created by the methods 
described above, a preliminary investigation of the candi- 
date flutter testing techniques was made. The analysis presented 
in Sections 6.3.1 through 6.3.3 is based on no external disturbance, 
that is, the only input is through the control surface 
excitation. Responses with turbulence as a noise source 
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FIGURE 6.8 FLUTTER TESTING CONTROL SURFACE INPUTS 
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are examined in Section 6.3.4 in order to evaluate the per- 
formance of the FSS in a noise contaminated environment. 

One concept in particular which was evaluated was the 
determination of the unaugmented characteristics of the wing from 
the closed loop response. Using the clean (noise free) data, 
the ideal circumstances necessary to test the feasibility of 
the concept are created. The interest in this concept stems 
from the fact that if it proves feasible, the open loop 
flutter boundary can be determined accurately with the FSS 
operating. The danger associated with flying with the system 
off while close to the unaugmented boundary could be avoided. 

An unpublished spectral analysis computer program, de- 
veloped under NASA DFRC sponsorship, was modified to include 
all the flutter methods mentioned in Section 6.1. This pro- 
gram, called SPA, created the plots for the analysis used 
in the next sections. 

6.3.1 Subcritical Open-Loop . - The subcritical unaug- 
mented simulated accelerometer response at wing station WBL 79 
(rear spar, 2.01 meters from the fuselage centerline) is 
shown in figure 6.9. The flight condition is M = 0.9 and 
h = 6096 meters. The control surface movement and command 
are also shown. The gust input is zero. Individual modal 
responses from this case and subsequent cases are presented 
in Appendix C.2. The increased response as the surface 
sweeps through the first wing bending frequency is evident. 

Figure 6.10 presents the results of the various flutter 
testing methods. The exact values are the eigenvalues of the 
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system of equations representing DAST I. The PSD of the 

response due to the doublet is seen in figure 6.11, The 

first and second mode (q and q ) are well defined. Damping for 

1 2 

mode q cannot be estimated because of the lack of a well defined 
3 

frequency peak. The estimate for damping in the second wing 

mode is accomplished by projecting the peak shape downward 

to estimate the half-power frequencies . (reference 12) 

Figure 6.12 is a plot of the real and imaginary parts 

of the H transfer function. Only the 10-30 hertz frequency 
2 

range is plotted because the lack of excitation outside this 
range creates erroneous or very little response in the low 
and high frequencies. Much more information is evident in 
the expanded frequency scale so better frequency and damping 
estimates can be extracted. 

The Amplitude-Phase plots shown in figure 6.13 indicate 
the three modes near the flutter frequency very distinctly. 

The square of H is shown in figure 6.14, The third mode 
is much better defined in this plot than in the PSD plot 
of figure 6,11. Too little information is present in the 
19.2 hertz mode to make a proper damping estimation. 

The inverse Fourier transform of H results in the im- 

2 

pulse response function shown in figure 6.15. The determi- 
nation of the damping estimate from the logarithmic decrement 
of the consecutive peaks assumes only one mode is present in 
the response. Recall from the previous plots that the 
first wing bending mode has much more response evident than the 
other two modes so it should dominate the initial part of the 
impulse response. This creates in essence the appearance of 
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the mode as the only mode present in the initial response. 

After approximately 0.3 seconds, the impulse response represents 
essentially useless information. 

6.3.2 Subcritical Closed-Loop 

The subcritical augmented response is shown in figure 6.16. 
The flight condition is the same as in the open loop case (h=6096 
meters and M=0.9). The affect of comparing this augmented 
response with the open loop response is the same as turning the 
FSS on and off in flight. The closed loop operation obviously 
has a large affect on the surface deflection during the doublet 
and sweep commands as can be seen by comparing figures 6.9 and 
6.16. 

The results of the application of the flutter testing 
methods on the augmented subcritical response is summarized 
in figure 6.17. Again three modes were identifiable but one 
is different than in the open loop case. The fuselage mode, 
q 2 , at 19 hertz did not appear but the third wing mode, q 4 , 
at 32 hertz was evident. 

The closed loop input spectrum of the surface movement due 
to the command sweep is shown in figure 6.18. The control 
surface is responding to the excitation of the first three wing 
modes. The first wing bending mode, q_^, has been moved to a 
lower frequency than was evident in the open loop case. The mode 
q^, the first wing torsion, is basically at the same frequency, 
however, it now is the mode that dominates the power spectrum. 
This implies that much more response should be evident in the 
mode q than at either the q or q wing mode frequencies. 
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The fuselage mode, I2' totally suppressed. 

The PSD plot of the wing response to a command doublet 
is shown in figure 6.19. The figure indicates the 8th struc- 
tural mode at 77 hertz besides the first three wing modes. 

The fuselage mode is not excited by the commands in these 
subcritical closed loop tests. Although the command is for 
a 25-hertz doublet, the PSD of the surface movement in figure 
6.20 does not show the same characteristics as found for the 
open loop case of figure 6.11. It is simply responding to 
the wing motions. 

Figure 6.21 is a Co-Quad plot of the closed loop transfer 
function. The estimate of the frequency and damping of 
was made from this plot, whereas the lower frequency modes 
were estimated from figure 6.22 which is an expanded scale 
plot of figure 6.21. No modes above q^ appear. 

The Amplitude-Phase plots of figure 6.23 were used to 

estimate the characteristic of q while the lower modes were 

4 

estimated from figure 6.24. The "flatness" of the phase 
angle plot makes estimates of q difficult. This flattening 
effect has appeared in wind tunnel testing of FSS models. It 
is an indication of a highly damped system. 

The modulus of the transfer function did indicate the 
presence of the q mode, but the peak was not usable as an 
estimator for damping. Figure 6.25 shows an expanded scale 
plot of | H 2 1 . The first wing torsion mode, q^ , is very well 
represented and indicates much more response in the third 
mode than the first. 
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As indicated from |H*| we would expect to see the 
mode dominate the initial part of the impulse response. 

Figure 6.26 verifies this fact. This response is not well 

enough defined to estimate damping. 

It has been seen that the FSS subcritical affects at 
6096 meters is to reduce the frequency of the wing bending 
mode, q , and increase its damping ratio. The fuselage mode 
appears totally suppressed and the first wing torsion mode, 

q , has had its damping reduced, although it is still suffi- 

3 

cient. There is some coupling of the FSS system to the 
fourth mode as evidenced by an increased response to the 
commanded sweep. 

The same testing methods were applied to try to deter- 
mine the open loop vehicle characteristics from the closed 
loop response. Surface movement is used as the input instead 
of surface commands. The results are summarized in figure 
6.27. No PSD plot is available because with the PSD method 
only the response is evaluated. The input is assumed to be 
white noise. Figure 6.28 shows the Co-Quad plots of the 
open loop transfer function. The three modes which appeared 
in Section 6.3.1 are evident. The q j mode peak is slightly 
noisy but a single peak is easily assumed. Comparing this 
plot to figure 6.12 demonstrates that the open loop charac- 
teristics have been reproduced. It should be noted that the 
q mode was indicated much better by determining the open loop 
characteristics from the closed loop response than from the 
basic open loop response. This was evident in both the 
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Co-Quad and Amplitude Phase methods. It is best illustrated 

by examining the broad spectrum plots of Amplitude-Phase 

of the open loop from the closed loop response seen in 

figure 6.29a and the open loop analysis seen in figure 6.29b. 

The 32 hertz mode, q , is represented much better in figure 

4 

6.29a than in figure 6.29b. 

The Amplitude-Phase plots for the open loop identifi- 
cation from the closed loop response is shown in figure 6.30 
for an expanded scale. The fuselage mode is faintly evident 
from this plot. 

The transfer function modulus shown in figure 6.31 

only indicates q and q . Only the first mode is sufficiently 

1 3 

indicated to estimate damping. Figure 6.32 is the impulse 
response function of the open loop system determined from the 
closed loop response. This, as much as any other plot, indi- 
cates that the subcritical open loop system is represented 
in a closed loop response. Comparing figure 6.32 to the 
impulse response function generated from the open loop re- 
sponse in figure 6.15 shows that the initial response is 
duplicated. 

6.3.3 Supercritical Closed-Loop . - The supercritical 
flutter test region conditions for this section are M = 0.9 
at h = 3048 meters. The time history of the simulated re- 
sponse is seen in figure 6.33. The first wing torsion mode 

q is the one which flutters first in this region. The time 

3 

history traces show that the accelerometer traces do not in- 
dicate as much response in q^ as it does in q^ during the 
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frequency sweep. 

The estimates of the closed loop frequencies and damping 
by the various methods are shown in figure 6.34. Figure 6,35 
shows the altered spectrum of surface movement to command 
input for the frequency sweep. The wing mode at approxi- 
mately 32 hertz is the flutter mode but its damping ratio 
with the FSS operating based on analysis is .92. So this 
mode was never present in the data. Only the q^ and q^ 
modes are evident as can be seen by the PSD plot of the re- 
sponse to a command doublet in figure 6.36. The 18.5 hertz 
q wing mode and the fuselage mode q^ are evident with a 
slight indication of the 77 hertz wing mode, q g , appearing. 

The PSD of the surface movement to the doublet command is 
almost identical to the response PSD. 

Both the Co-Quad and Amplitude-Phase plots shown in 
figures 6.37 and 6.38 respectively, are almost representa- 
tive of classical examples of the methods. For the two 
modes excited, the information for the frequency and damping 
estimates is well defined. 

The transfer function modulus seen in figure 6.39 
would indicate that the impulse response function should be 
dominated by the 21 hertz fuselage mode q^. Figure 6.40 shows 
this to be true although a damping estimate cannot be made 
from the impulse response. 

The FSS completely suppresses the known flutter mode 
q and good damping is indicated in the wing bending mode 
q . The response in the fuselage mode has increased sub- 
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stantially however. 

The attempt at extracting the open loop characteristics 
of the wing from the closed loop supercritical response is 
summarized in figure 6.41. The open loop roots for the 
first three modes range from 19.3 to 19.7 hertz. The Co-Quad 
plots in figure 6.42 indicate a lightly damped mode at about 
19.6 hertz. This is probably the fuselage mode q 2 since 
is highly damped and q^ is unstable for the open loop case. 

The source of the peaks below 14 hertz are unknown since 
very little input excitation exists at these frequencies 
it is probably just noise due to an extremely small denom- 
inator term. The Amplitude -Phase plot of figure 6.43 
indicates very little information. The Phase plot does 
indicate the mode but little else. The transfer function 
modulus presents inconclusive results as is evident in 
figure 6.44. From figure 6.44 one would expect the impulse 
response function, shown in figure 6.45, also to be incon- 
clusive. The plot is representative of a highly damp system 
of indeterminate frequency. 

Although some of the open loop modes are present in 
the supercritical response, the modes involved in flutter 
cannot be reliably identified, 

6.3.4 FSS Performance In Turbulence .- The preliminary 
analysis on flutter testing methods performed in the 
previous section was done on a noise-free data. This allowed 
the evaluation of the concept of determining open loop 
flutter characteristics from the closed loop responses. 

In a flight test environment noise sources from turbulence. 
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engine noise, and separated flow will cause contamination 
of the data. The methods described previously all determined 
the frequency response transfer function in a way which will 
eliminate noise in the data. This is only true if the 
response due to the command excitation is a large portion 
of the total response. If on the other hand the majority 
of the response is due to the noise inputs the determination 
of the transfer function becomes an extremely difficult 
task (if not impossible) . 

Figure 6.46 shows the transient response of the ARW-1 
FSS at M = 0.9, h = 3048 meters, and with a .3048 meters/sec 
rms gust level in the vertical axis. Peak accelerations are 
approximately 6 g's. The peak surface deflections are 
approximately 5.5 degrees. Since the amplitude of the 
command inputs are only 1 degree, it should be expected that 
with this input magnitude a reasonable estimate would be 
difficult to achieve. The "bursting" in the accelerometer 
trace is considered representative of the traces from wind 
tunnel and flight test results. 

A determination of the rms level of the surface de- 
flection^ , in figure 6.46 has not been completed to date, 
s 

Therefore no comparison in surface deflection requirements 
to BWC ' s predictions were made. 
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CHAPTER 7 


CONCLUDING REMARKS AND RECOMMENDATIONS 

7 . 1 Concluding Remarks 

7.1.1 Vehicle Development 

1. Based on analysis, the Boeing FSS design is 
predicted to stabilize the ARW-1 flutter mode 
the required 20 percent margin past the open loop 
boundary. 

2. A digital remotely augmented rigid body control 
system which provides attitude and load limiting 
protection was designed and evaluated for the 
DAST I. Based on wind tunnel and DATCOM es- 
timated rigid body aerodynamic characteristics, 
this RAV system meets the design specifications 
and has performed satisfactorily during real- 
time piloted simulations. 

3. An on-board backup autopilot control system was 
designed to provide an envelop limiting command 
capability from an F-104 chase aircraft. The 
design has flight path control and wing leveling 
as automatic functions. Linear analysis and 
real-time simulations show this system to be 
satisfactory. 

7.1.2 Evaluation of Flutter Methods 

Simulated data was generated from a state variable 
representation of the DAST ARW-1 equations of motion. 
This data was used to evaluate candidate flutter 
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testing methods. The following conclusions can 
be drawn from the evaluation: 

1. The five candidate methods all predicted modal 
frequencies accurately. 

2 . The damping estimates for noise free 
data are good for all the methods. How- 
ever the PSD and transfer function modulus 
are more optimistic than the others. 

3. If one mode dominates the accelerometer response 
the impulse response function can 

identify that one mode well. 

4 . Open loop and closed loop characteristics are both 
identifiable from their respective responses. 

5. The concept of determining open loop 
characteristics from the closed loop responses 
with the flutter testing methods works very well 
for subcritical data. Although some modes are 
present in determining the open loop character- 
istic from the supercritical responses using 
Co-Quad and Amplitude-Phase methods, the open loop 
flutter mode was not readily identifiable with 
any of the methods examined. 

7 . 2 Recommendations 

1. Further investigation of determining unaugmented 
modal characteristics from closed loop responses 
should be investigated. The degradiation in the 
frequency and damping estimates by this method with 
noise input needs to be examined. 
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2 . 


The effect noise has on the five techniques 
examined for near real-time testing needs exam- 
ination. The results obtained from "noise free" 
data could be quite different if the same study 
were conducted with noisy inputs. 

3. The non-linear loaded transfer function of the 

FSS actuator should be included in the simulation 
to make an accurate assessment of the surface 
rms displacement and rate requirements in turbulence 
and check for possible limit cycle tendencies. 
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APPENDIX A 


PAST ARW-1 AERODYNAMIC CHARACTERISTICS 

This appendix presents the rigid body aerodynamic data 
developed for the DAST I flight simulator and control systems 
analysis. The derivatives for the DAST are given as func- 
tions of Mach number and angle of attack to make them com- 
putable with the NASA DFRC simulation. The lateral-direction- 
al derivatives are in the body axis coordinate system while 
the longitudinal derivatives, C L , C D , and C m are in the sta- 
bility axis system. 

A. 1 Static Aerodynamic Characteristics 

A. 1.1 Longitudinal Coefficients - Lift, drag and pitching 
moment coefficients for the Mach number conditions tested 
in the NASA LaRC 8-foot Transonic Dynamic Tunnel are shown 
in figure A.l. The angle of attack of the wind tunnel tests 
varied from approximately -4 degrees to 12 degrees. The 
pitch-up characteristics of the DAST I is evident in these 
figures. When M = .90, the change in lift curve slope and 
pitching moment coefficient becomes very pronounced at about 
seven degrees angle of attack. This behavior was present on 
the F-8 SCW aircraft and is attributed to spanwise separation 
of airflow near the wing tip. 

In order to obtain the angle of attack derivatives 
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PITCHING MOMENT COEFFICIENT 
ANGLE OF ATTACK (PAGE 1 OF 6) 










C D and C m in the next section, trim points from the C L 

a a 

versus C curves in figure A. 2 were found. This allows 
m 

a determination of the trim C L and C D from the drag polars 
also presented in figure A. 2, referring these conditions 
back to figure A. 1. The slope of the drag and pitching 
moment curves can then be established at trim conditions. 


A. 1.2 Angle of Attack Derivatives - The plot of C L 

a 

at various Mach numbers and angles of attack is shown in 

figure A. 3 along with the effective intercept coefficient 

C T . The maximum value of C T is 7.74 per radian (.135 per 
° a 

degree) and occurs at M=.98 and two degrees angle of attack. 
The dashed lines after M=.98 represent the extrapolated data. 
By using an effective slope intercept coefficient C L ^, non- 
linear lift curves can be represented. Observing C L for six 

a 

degrees and eight degrees angle of attack in figure A. 3, the 
large change in C L associated with the tip separation can be 


The power-off rate of change of drag coefficient with 

angle of attack is shown in figure A. 4. Here C Dq is also 

the effective intercept coefficient. These derivatives 

are determined at trim conditions. It is interesting to 

note that C is still negative at almost two degrees angle 
D a 

of attack. A great amount of testing will be done below 

two degrees angle of attack in the FSS flight experiment. 

The dip in C D and C D corresponds to the dip in C L 
a ° a 
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FIGURE A. 2 STATIC STABILITY DATA AND FIREBEE DAST DRAG POLARS 
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FIGURE A. 2 (CONTINUED) (PAGE 3 OF 6) 
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FIGURE A. 2 (CONTINUED) (PAGE 6 OF 6) 
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and Cl q associated with the tip separation problem dis- 
cussed earlier. 

The rate of change of pitching moment coefficient 
with angle of attack is shown in figure A. 5. The large 
positive values of C ma due to pitch-up are very apparent. 
It can be seen that the change in sign of C ma occurs very 
abruptly. The slope intercept coefficient Cn^ is also 
shown in figure A. 5. 


A. 1.3 Longitudinal Control Derivatives - Values of 

the rate of change in lift coefficient with horizontal 
tail deflection are shown in figure A. 6. Control surface 
effectiveness increases substantially between M = .9 


and M = .98. In comparison with the basic Firebee drone, 
Cl 6e for DAST is approximately 27 percent greater although 
the ratio of Firebee reference areas to DAST reference 


area is 1.07. 

The rate of change of drag coefficient with hori- 
zontal tail deflection is presented in figure A. 7. Trim 
points were found from the Cl versus C m curves of figure 
A. 2 and the increase in drag coefficient at that hori- 
zontal tail deflection is determined from the correspond- 
ing drag polars also in figure A. 2 The increment in drag 
coefficient is then divided by the amount of tail deflec- 


tion to trim the aircraft and the result is Cp. . A 

°E 

"hump" in the C D curve in figure A. 7 is evident around 
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FIGURE A. 5 VARIATION IN RATE OF CHANGE OF PITCHING MOMENT 

WITH MACH NUMBER AT VARIOUS ANGLES OF ATTACK 
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FIGURE A. 5 (CONCLUDED) 
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FIGURE A. 7 VARIATION IN RATE OF CHANGE OF DRAG COEFFICIENT WITH 

HORIZONTAL TAIL DEFLECTION WITH MACH NUMBER AT VARIOUS 
ANGLES OF ATTACK FOR THE DAST 



is attributed to a 


M = .80. The initial rise in C D ^ 

E 

rearward shift in the aircraft aerodynamic center . The 
decrease is associated with less trim deflection requir- 
ed because of the separation at the wing tip. 

The derivative C m is presented in figure A. 8. All 

III p 

° E 

wind tunnel data were referenced to 25 percent MAC of the 
wing. The magnitude of C m6 for DAST is approximately 
twice that of the basic Firebee drone. 

A. 1.4 Lateral-Directional Coefficients - Wind tun- 

nel data for the DAST I showed that the lateral-direction- 
al coefficients, C^, C n , and Cy were mildly non-linear with 
respect to sideslip angle. No provisions exist in the sim- 
ulation equations to include this effect. Since the effect wors 
ened with larger sideslip angles the assumption was made 
to use the values which exist at small (+3°) sideslip. The 
validity of this assumption rests in the fact that the 
sideslip angle, 8, will be nearly zero in flight. Large 
rudder inputs are required to increase it substantially. 

The variations in C^, C^, and Cy with angle of at- 
tack and Mach number is shown in figure A. 9 . Only at 
M = . 95 does there appear any change in Cy with angle 
of attack. Since C y versus angle of attack at 8 = 2° 
and M = .95 does not show this effect, it is assumed 
that Cy was independent of angle of attack. However, 
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FIGURE A. 8 
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VARIATION IN RATE OF CHANGE OF PITCHING MOMENT COEFFICIENT WITH 
HORIZONTAL TAIL DEFLECTION AT VARIOUS MACH NUMBERS AND ANGLES OF ATTACK 
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FIGURE A. 9 VARIATION IN ROLLING MOMENT, YAWING MOMENT, AND 
SIDEFORCE COEFFICIENTS WITH ANGLE OF ATTACK AT 
VARIOUS MACH NUMBERS AND SIDESLIP ANGLES 

(PAGE 1 OF 6) 
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C-£ and C n do show a dependency on angle of attack at all 
Mach numbers. Rolling moment coefficient, C £ , seems par- 
ticularly affected in the region where the pitch-up 
characteristic is observed. 

A. 1.5 Sideslip Derivatives - The rate of change in 

side force coefficient with sideslip angle, Cy^, is pre- 
sented in figure A. 10. The vertical tail component is 
shown because of its importance in estimating some of 
the dynamic derivatives in Section A . 2 . 2 . The side 
force produced by the vertical tail is assumed to be the 
same on both the DAST and basic Firebee vehicles. The 
differences in C v seen in figure A. 10 is due to the 
different reference areas of the two aircraft. 

Figure A. 11 presents the directional stability deriv- 
ative, C„ . Values for the DAST are approximately 60 

n e 

percent of the Firebee' s magnitude. The differences 
are attributed to the wing-body interference affects at 
low angles of attack and wing tip separation at high 
angles of attack. 

The rate of change of the rolling moment coefficient, 
Ci^r due to sideslip angle is shown in figure A. 12. 

The large jump in C£ g which occurs around eight degrees 
angle of attack above M = .80 is caused by the increased 
spanwise separation of the downstream wing tip. This 
results in an effective dihedral of very large proportions. 
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MACH NUMBER 

FIGURE A. 10 VARIATION IN RATE OF CHANGE IN SIDEFORCE COEFFICIENT DUE TO 
SIDESLIP ANGLE WITH MACH NUMBER 
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FIGURE A. 11 VARIATION IN RATE OF CHANGE OF YAWING MOMENT COEFFICIENT DUE TO 
SIDESLIP ANGLE WITH MACH NUMBER AT VARIOUS ANGLES OF ATTACK 
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MACH NUMBER 

FIGURE A. 12 VARIATION IN RATE OF CHANGE OF ROLLING MOMENT COEFFICIENT 
DUE TO SIDESLIP ANGLE WITH MACH NUMBER AT VARIOUS ANGLES 
OF ATTACK 




The low angle of attack values of C t for DAST are close 
to the values of the Firebee. The magnitude of at 

low angles of attack are also close to those experienced 
on lifting bodies. 


A. 1.6 Lateral-Directional Control Derivatives - The 

rate of change in sideforce coefficient with rudder de- 
flection is presented in figure A. 13. These values were 
calculated from wind tunnel data using rudder deflections 
from zero degrees to -IQ degrees. From this data it was 

observed that C v did not vary with angle of attack and 

y<S R 

was close to the Firebee values as one would expect. 

Figure A. 14 shows the variation in yawing moment 

coefficient due to rudder deflection. The values for 

Cn,. in this figure are approximately 60 percent of the 
o R 

Firebee' s values. This change can be attributed to the 
large reference span of ARW-1. 

The variation of rolling moment due to rudder de- 
flection, Ci , is shown in figure A. 15. It is appar- 
6 R 

ent that at large angles of attack the value of Co is 

OR 

reduced substantially. Otherwise, the magnitude is re- 
presentative of those experienced on the basic Firebee 
drone. 


The rate of change in sideforce coefficient due to 

differential tail deflection, C v , is presented in fig- 

Y $d 

ure A. 16. The magnitude of C Vr for DAST is five times 

6 d 

that of the standard drone. This large difference is 
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FIGURE A. 13 VARIATION IN SIDEFORCE COEFFICIENT DUE TO RUDDER 
DEFLECTION WITH MACH NUMBER 
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FIGURE A. 15 VARIATION IN ROLLING MOMENT COEFFICIENT DUE TO RUDDER DEFLECTION 
WITH MACH NUMBER AT VARIOUS ANGLES OF ATTACK 





considered to be due to the different downwash effects be- 
hind the thin uncambered supersonic wing of the Firebee and 
the much thicker ARW-1. Differences in reference span also 
contribute. 

Figure A. 17 shows the variation in yawing moment co- 
efficient due to differential tail deflection, C n , • The 

a 

Firebee values for C n ^ are approximately sixty percent of 
the DAST values at higher Mach numbers. Flight test results 
of the basic Firebee II drone indicate values of the same 
magnitude but of opposite sign of the wind tunnel values 
compared here. 

The rate of change in rolling moment coefficient due 

to differential tail deflection, C ff . , is shown in figure 

6 d 

A. 18. The high angle of attack region appears to influence 

this control derivative. The wind tunnel data indicated 

that the variations in C 0 due to collective differential 

^d 

tail deflection and elevator deflection are insignificant. 

The magnitude of C» on DAST is about forty percent higher 

°d 

at low angles of attack than that experienced with the 
Firebee. 

A. 2 Dynamic Derivative Estimates 

All derivatives in this section were determined by the 
latest USAF DATCOM (reference 6) methods and those of ref- 
erence 13. Static wind tunnel data were used as input in- 
formation to the empirical expressions whenever appropriate. 
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FIGURE A. 17 VARIATION IN YAWING MOMENT COEFFICIENT DUE TO DIFFERENTIAL 
HORIZONTAL TAIL DEFLECTION WITH MACH NUMBER AT VARIOUS 
ANGLES OF ATTACK 


//// 






*T COEFFICIENT DUE TO DIFFERENTIAL HORIZONTAL TAIL 
SR AT VARIOUS ANGLES OF ATTACK 


o o 




A. 2.1 Longitudinal Dynamic Characteristics - The 

estimated values for lift coefficient due to rate of 

change of angle of attack are presented in figure A. 19. 

It is assumed that the horizontal tail increment was the 

only significant contribution to the derivative. For 

the horizontal tail contribution, C T . , is calculated by 

■““H 

the expression: 


: ^ H " 2 CL “H ^ ^ £ 


(A. 1) 


The lift curve slope for the horizontal tail was obtain- 
ed from reference 3, page 29. The term V R is the hori- 
zontal tail volume coefficient defined as: 



S H 

S 


(A. 2) 


Where X„ is the distance from the airplane eg to the 
horizontal tail aerodynamic center, S H and S are the tail 
and wing areas respectively, and c is the mean aerody- 
namic chord of the wing. 

The rate of change of downwash angle with angle of 
attack, de/da, is shown in figure A. 20. It was deter- 
mined by the following equations: 


de 

da 


M 



M=o 




{m 


lM=o 


(A. 3) 
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de 

dot M _ 0 


= 4.44 


1+AR 1 * 7 




h H^ 


aI 2 *h 






cos A 


C/4 


(A, 4) 


where 

A = taper ratio 

AR = aspect ratio 
b = wing span 
A C /4 = < 3 uarter chord angle 

«, H = distance from 25% c of the wind to 25% c of the 
horizontal tail 

h H = vertical distance of the horizontal tail from 
the root chord plane of the wing 
The ratio of lift curve slopes used in equation (A. 3) 
was determined from the wind tunnel data. Since the 
transonic force-break Mach number for the DAST I was ap- 
proximately M = .9 (reference 6, figure 4.1,3.2~53b) f the 
wing-body lift curve slope for M = .40 was used in the de- 
nominator of equation (A. 3). 

The dynamic pressure ratio at the horizontal tail, 
n H r used in equation (A.l) is shown in figure A. 21. 

This ratio was estimated by s 



(A. 5) 
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FIGURE A. 21 VARIATION IN HORIZONTAL TAIL DYNAMIC PRESSURE RATIO WITH MACH NUMBER 
AT VARIOUS ANGLES OF ATTACK 


o o 



(A. 6) 



c 


Z — Xjj tan ( e - a ) 

(A. 8) 

1.62 C L 

e — 

(A. 9) 

n ar 


z = 0.68 c aJ c + • 15 ] 

(A. 10) 

W V Dq V c ) 



The term C n in equations (A. 7) and (A. 10) does re- 

u o 

present the conventional drag coefficient at zero angle 
of attack, and C L is the total lift coefficient at the 
particular flight condition of interest. Both Cj-, o and 
C T were determined from the wind tunnel data. 

The values for C-^ in Figure A. 19 are not signifi- 
cantly different than that of the Firebee at low angles 
of attack but degrade at high Mach numbers at high angles 
of attack. 

The derivative C L ^ is composed of a wing contribu- 
tion and horizontal tail contribution. The wing contri- 
bution is estimated by: 

c r 

(A. 11) 


AR + 2 cos ^c/4 
AR B + 2 cos A c /4 
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B 


(A. 12) 



M 2 COS 2 

fl/2 + 


A c/i+ 



M=o 


(A. 13) 


where X w is the distance from the airplane eg to the wing 
aerodynamic center. 

The tail contribution was found from : 



2 C T 


n H 


V, 


‘HI 


H 


M 


(A. 14) 


Summing equations (A. 11) and (A. 14) yields the values 
plotted in Figure A. 22. Values of C L for the DAST and 

T. 

Firebee are comparable. 

No explicit methods exist to calculate the wing 

contribution to the derivative C m . except for triangular 

a 

wings. For most conventional airplane configurations 
its contribution is small. For the horizontal tail con- 
tribution the estimate is: 


- -2 c L n H v H %* (A.i: 

a a H U H c da 

The derivative C m> is plotted in Figure A. 23 as a 

a 

function of Mach number and angle of attack. At low 

angles of attack C m# for DAST I is about twice that of 

a 

the basic vehicle. Above M = .8 and at high angles of 

attack C decreases sharply. Improved methods were 
in * 

a 

used in the DAST calculation. 

The variation in pitching moment coefficient with 
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FIGURE A. 23 VARIATION IN PITCHING MOMENT COEFFICIENT WITH RATE OF CHANGE 
ANGLE OF ATTACK 


change in pitch rate is considered the sum of the wing 
and tail contributions. The wing contribution is found 
by : 






AR 3 tan 2 

A c/4 j. 3 

AR*B+6 cos Ac/4 B 

AR 3 tan 2 

^c/4 + 3 

AR+6 cos 

A c/4 


(A. 16) 


'”' qwb lM=o 


= .71 C 


a, 


cos A 


*w 


c /4 


AR 


[2 (hi V + 1/2(M 
l \ C ) \ C ) 


AR + 2 cos A c /4 


24 


AR 3 tan 2 Ac/4 + 1 
AR + 6 cos A c /4 8 

(A. 17) 


In equation (A. 17), C £ is the spanwise average value 

_ a w 

of the wing section lift curve slope. 

The tail contribution to C m is estimated from: 



■ ' 2 Cb « H " H V H 



c 


(A. 18) 


The derivative C m is presented in Figure A. 24. 

ra q 

The DAST values for Cxnq are about 13 per cent higher than 
those of the Firebee at the same angle of attack. The 
increased aspect ratio for the DAST works to its bene- 
fit. 


A. 2. 2 Lateral-Directional Dynamic Characteristics - 
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The sideforce coefficient due to angular roll velocity 
is plotted in figure A. 25. This derivative is composed 
of two components, one for the vertical tail and one for 
the wing. At low angles of attack the vertical tail is 
the main contributor while at high angles of attack the 
wing component is the larger of the two. 

The vertical tail component is estimated from: 


Ct, = 2 Zv cos 01 " *v sin a C v _ (A. 19) 

y P v b 

where Z v is the distance from the vertical tail aerody- 
namic center to the body x-axis and £ v is the distance 
from the vertical tail aerodynamic center to the aircraft 
center of gravity. 

The wing-body component is estimated from: 

Cy Pwb] [ Cy Pwb ] . r AR + 4 COS A c/4 1 j- AR-B + cos A c /4 ~ | 

C L / M \ C L / M=0 lAR.B + 4 cos A C /4J l AR + cos A c/4 J 

(A. 20) 


the value of ( C v / C L ) 
y Pwb 

3a of reference 6. 


M=o 


is found from figure 


7. 1.2.1- 


The variation in yawing moment coefficient with roll 
rate is estimated from its wing and vertical tail compon- 
ents. For a wing with no twist or flaps its contribution 
is estimated from : 




u 

M 



(A. 21) 
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where : 


0Cj 

IT 


!J2. K/B 


(A. 22) 


and ( BC Jp /K) is found from figure 8«1 of reference I 3 
The term K is the ratio of average wing section lift 
curve slope, C» , to 2n, and 3 =V 1-M 2 . 


The roll damping derivative, C^, in equation (A. 21) 

is determined later in this section. The term (Cn p /CL) Ct,=o 

M 

was determined by : 




Cl=o 

M 


A R + 4 cos A c/ / 4 


AR’B + 4 cos A c /4 
C 


AR-B + % (AR-B + cosAc/ 4 ) tan 2 A c /4 
AR + * 5 (AR+ cos A c / 4 ) tan 2 A c / 4 


'n r 


C L 


C L =o 

M=o 


(A. 23) 


where 


3e\ 


[jIE 

' C - L. 


C L =o 

M=0 


1/6 


. fx tan A c / 4 tan 2 I^c/^ 
AR+ 6 (A R + cosA c/4 ) 1 1 Al C -^- + n ) 


A R + 4 cos A c / 4 


CA.24) 


The vertical tail contribution can be found from 

r Z v cosa - sina^ 


C npv = - 2/b U v cosa + Z v sina)^- 


. V — j C Y 
b / "v 


(A. 25) 
3. 


The summation of equations (A. 21) and (A. 25) is present- 
ed in figure A. 26. 

This estimated derivative differs significantly 
from the basic drone. This variance is attributed to the 


298 




.80 1.0 1.2 
H NUMBER 

IFFICIENT DUE TO ROLL RATE WITH 
OF ATTACK 



very different wing geometries of the two vehicles. 

The roll damping derivative is derived from contribu-v 
tions of wing, horizontal tail, and vertical tail. The 
wing contribution at zero lift was estimated from equation 
(A. 22). When C L is greater than zero, the additional wing 
component is estimated by: 

C L 2 1 1 + 2 sin 2 Ac /4 


AC* 


Pw 


1 

8 


n AR cos 2 Ac/4 


AR + 2 cos A 


c/4 


AR + 4 COS A c / 4 




;p 
. 1 

1 + ^Do 




(A. 26) 


where C n is the drag coefficient at zero angle of attack. 
u o 

The horizontal tail component is approximated by: 

12 


= .5 


'PH 


(c. 

y s h 

^e\ 

V* , i 

V b REF 

\ b l 


(A. 27) 


where C» is the tail contribution based on its own geo- 

V *P / h 

metry and found from equation (A. 22) . 

The vertical tail contribution is estimated from: 


'Pv 


= 2 


J v 
"E - 'i 


J 


'*8 


(A. 28) 


v 


Figure A. 27 shows for the total aircraft. The high- 
er aspect ratio of DAST I causes to be about twice that ob- 


served on the Firebee. 

The sideforce coefficient due to angular yaw velocity is 
estimated from a vertical tail component only. The wing con- 
tribution is considered negligible. This derivative may be 
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found from: 


'Yr 


C 


Yr 


v 


= - 1 ( 
b V 


«, v cosa + Z v sina j C 


Ye 


(A. 29) 


v 


It is shown in figure A. 28. The much larger reference 
span of ARW-1 causes Cy r to be only sixty percent of the Fire- 

bee values. 

The yaw damping derivative, C nr , is estimated from con- 
sideration of two components, that due to the vertical tail 

and that due to the wing . 

The wing contribution can be estimated by: 

(A - 30) 


- = 

' nr W \c^ j 


c L 2 + 


C Dt 


where the ratio (C nr /C L 2 ) is found in fi 9 ure 9 * 4 of reference 
13, and the ratio (C nr /C Do ) can be found from figure 9.5 of 
reference 13 . The term C Dq is the drag coefficient at zero 
angle of attack and zero control surface settings. 

The vertical tail component is found from: 


C nr ^ = -1 cosa + Z v sina j 2 Cy^ (A. 31) 

The larger values for C nr , seen in figure A. 29, for the DAST I 

vehicle are due to the larger span. 

The rolling moment coefficient due to yaw rate is present- 
ed in figure A. 30. This derivative can be estimated from a 

wing and vertical tail component. 

The wing component of the derivative for a zero-twist, 
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FIGURE A. 30 VARIATION IN ROLLING MOMENT COEFFICIENT DUS TO YAW RATE WITH 
MACH NUMBER AT VARIOUS ANGLES OF ATTACK 




un- 


■flapped, and no dihedral wing is estimated by: 
tn 

^ o 


c * r ~ C L 
A w 


(A. 32) 


c L 


C L=0 

M 


where the slope of the rolling moment due to yaw rate at zero 
lift is given by: 


Cl 


ar(1 . b 2) + AR + 2cos A c/4 . tan 2 Aq/Jl 

1 + 2B(AR*B + 2cos A c / 4 ) AR' B + 4cos A c /4 8 


Cl=0 

M 


+ AR + 2 cos A c / 4 
AR + 4 cos A c / 4 


tan 2 A c/4 


and 


(A. 33) 


I Cl=0 
m=0 


\ Cl /c l =o 

M=0 


xs 


obtained from figure 9.1 in reference 13 as 


a function of quarter chord sweep angle, taper ratxo, and 
aspect ratio. 

The contribution of vertical tail was found from: 


C* 


r v 


_ («, cosa + Z v sina) (Z v cosa - & v Sina) 

b 2 


v 


(A. 34) 
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APPENDIX B 


PROJECT RESPONSIBILITIES AND TIMETABLE 

This appendix will identify the project responsibilities 
of the three participating organizations in the DAST I pro- 
gram. A milestone chart indicating the relative state of pro- 
gress toward the first flight is also shown. 

B. 1 Project Responsibilities 

The NASA Langely Research Center is the lead center in 
the DAST program. The NASA Dryden Flight Research Center has 
the responsibility for flight operations. The Boeing Wichita 
Company has the responsibility for design and fabrication of 
the flutter suppression system. Within these guidelines, the 
specific project responsibilities are identified in Table B.l. 
During flight testing BWC will participate as a consultant, 
but the principal investigators from DFRC and LaRC will be 
responsible for the conduct of the flight experiments. 

B . 2 Project Timetable 

The DAST I project milestones are listed in Table B.2 
The program began in July, 1974, and will continue through 
late 1979. First flight is planned for March, 1979, and the 
flight activity rate is expected to be one flight per month 
for the six scheduled flights. A system check flight (mile- 
stone number 26) with a standard wing is planned for late 1978 


307 


TABLE B.l PROJECT RESPONSIBILITIES 


Tasks 

LaRC 

DFRC 

BWC 

Project Funding and Procurement 

~ 1 


■ 

Define Goals and Objectives 

X 

X 

■ 

Flutter Experiment Requirements 

X 

X 


Wing Design and Construction 

X 


■ 

FSS Design 




FSS Fabrication 



X 

Wing Instrumentation 

X 


■ 

Flutter Instrumentation (PCM) 

X 

X 


Wind Tunnel Tests 

X 


X 

DAST Avionic Systems 


X 


DAST Simulation 


X 


DAST Flight Controls 


X 

■ 

Flutter Testing Methods 

X 

X 


Flight Planning 


X 

■ 

Flight Operations 


X 


Flight Safety 


X 


Data Collection 


X 


Data Reduction 

X 

X 


Flight Testing 

X 

X 


Reporting 

X 

X 

X 
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TABLE B . 2 - DAST I PROJECT MILESTONES (PAGE 1 OF 2) 


MILESTONES 


1 ARW - 1 Design 


6 Fabricate ARW-1 


7 Install Instrumentation 


8 FSS Preliminary Design 


9 FSS PD PDR at LaRC 

10 FSS Final Design 

11 FSS PDR 


12 FSS CDR 

13 Install ARW-1 for GVT 


14 Install FSS at LaRC 


15 RPRV Facility Check 


1974 


2 

Rigid Body 1 

W-T Tests 


3 

Develop FLT 

Simulator 


4 

PDR at LaRC 



5 

FLT Control 

Law Design 
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TABLE B.2 - DAST I PROJECT MILESTONES (PAGE 2 OF 2) 



APPENDIX C 

ARW-1 FLUTTER SIMULATION .. 

This appendix describes the development of the state vari- 
able representation of the ARW-1 equations of motion used in 
the FSS analysis. Transient responses at the accelerometer 
location at WBL(79) were generated using a modified version 
of the CONTROL program (reference 8) . This system of equations 
includes inputs from the control surface and turbulence from 
a Dryden gust model. 


C. 1 State Variable Formulation 

The Laplace representation of the ARW-1 equations of mo- 
tion was given in equation (4.2) of section 4.1.1. To gener- 
ate time histories of this system using CONTROL the equations 
were transformed to the form: 



This transformation was accomplished in the following manner 

J x i}= J gi ft)} < c - 2 > 

M" {*$■{’! *4 <c-3> 

The augmented states which represent the lag in buildup 


of circulation are: 






(C . 5) 



S + 


U 0 b 2 



(C. 6) 


(C . 7) 


In time domain the equations (C.4) through (C.7) become: 






(C . 8) 



(C.9) 



(C.10) 



(C.ll) 


If the higher wing modes are truncated and seven modes 
are kept , the number of states generated is 42. This trun- 
cation causes no problems since the accuracy in the higher 
modes is dubious due to the 400 rad/sec limitation on the 


aerodynamic approximations and the higher modes were shown to 
have no effect on the flutter characteristics. The x-trans- 
lation and z-translation rigid body modes were also dropped 
to reduce computer core memory requirements. They also were 
shown to be irrelavant to the flutter problem. The rigid 
body pitch mode was kept. 

The wing's response due to gust inputs also requires 
augmented states to represent the lag in buildup of circulation. 
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These terms are: 


X es _ 

4 3 S + U ogi 

X = s 

44 s + U 0 g 2 

X = s 

45 S + U c g 3 

X = 

46 S + U c g 4 


W g (S) 

V s 1 

w g (s) 

Wg(S) 


(C . 12 ) 
(C . 13) 
(C . 14) 
(C . 15) 


or : 


X = -U ogi X 43 + W g 

*T 0 


(C.16) 


X = -U Q g X + w 
44 ° 2 44 y 


CC.17) 


X - = '^3 X 45 + ^ 


45 


(C . 18) 


X = -U D g X + W 
46 ° 4 46 y 


(C . 19) 


The Dry den form of the spectra for the turbulence vel 
ocity (MIL - F8785B ASG) for a vertical gust is: 


w. 


w 


( L w fi ) 

T [i + (l w b)‘ ] 2 


L 1 + 3 

(&) = a £ — ; rz — TT- 


CC. 20) 


In converting from spatial frequency, equation (C.20) 
becomes: 
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Q = 


w 

U, 


<f> w (w) = a 


L 

_ ,2 w 


1 + 3 Mfe) 


(C . 21) 


,g ' W ^ [l + (Lw g-j 2 ] 2 

The frequency, to, is related to the Laplace variable by co= s/j 
so equation (C.21) reduces to : 


w„ ( s) 

2 L 

- a 2 — - 

1 + 

3! 

^±1 

n (s) 

« "Uo 

1 +1 

fe s/j ] 

2 - 


(C.22) 


J w 


” °w nV 


w 1 - 31 U, 


Hfe-r ] 


I ' / \ 

0 ,r» fi + v? s ) 

W V^o ( 1 + ^ s 


w 





Wg (S) _ WgTs) 
nTsT n (s) 


(C. 23) 


From equation (C.23) the expression describing the ver- 
tical gust is: 

" h + YT L ' 

'w\/7u7 


w 


./sT ( I . */T nr 1 i 




(C.24) 


After equation (C.24) is placed in a standard transfer 
function form the state variable representation of the vertical 
becomes : 
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+ x 48 

+ X 48 


X 


4 8 


Uo 


CT W 



r) 


(C . 25) 
(C . 26) 
[C .21) 


X 47 + 


a w 



n (c . 28) 


The input variable , n t , represents Gaussian white noise 
of unity variance. 

Although the control surface is considered one of the 
ARW-1 modes, it was considered separately so that actuator 
dynamics could be included in the equations. 


X 49 = 6g 



X51 = X50 - UobjX^ 


'52 


■ X 


50 


“ U O b 2 X 52 


X53 = X 50 


Uo b 3X 5 3 


(C . 29) 

(C . 3 0 ) 
(C . 31) 

(C . 32) 

(C . 33) 


X54 = X 50 ~ U Q b 4 X 5 4 


(C . 34) 
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Equations (C.31) through (C.34) represent the lag in 
buildup of circulation due to control surface movement. 
The transfer function representing the actuator is : 


6 s _ (846) 2 (C.35) 

6 c s 2 + 758s + (846) 2 

This is a simplification of the closed loop no-load trans- 
fer function described in the first critical design review 
at BWC held in October of 1977. The complete transfer func- 
tion is presented in Section 5.3. The poles at 10.14 radians 
and 375.4 radians were neglected because of pole-zero cancel- 
lation. The pole at 10.14 radians is cancelled by a zero in 
the actuator transfer function and the pole at 375.4 radians 
is essentially cancelled by the pole in the compensation 
filter shown in figure 4.5. The higher frequency second order 
poles in the complete actuator transfer function are neglected 
due to their frequencies being much higher than the frequency 
range of the wing modes . 

The equations of motion now become: (C.36) 
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Moving all the terms in equation (C.36) to the right 
hand side and multiplying by the inverse of the mass matrix 
( [M] “ 1 ) the vector is defined. This only leaves X; 0 

undefined in order to use equation (C.l) . 

Acceleration of the control surface is 



+ 6 (846) 2 

c 


Equation (C.37) differs from the control surface terms in 
equation (C.36) in that they represent the coupling of the 
control surface to the wing modes in equation (C.36) while in 
equation (C.37) the coupling of the wing modes to the control 
surface are represented. The coupling of the wing modes to 
the control surface were necessary to prevent a destabiliza- 
tion of the sixth wing mode in closed loop FSS operation. 

The inputs into the system are illustrated in figure C.l. 

The noise input to the gust is uncorrelated Guassian white noise 
of unity variance. 
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FIGURE C.l DAST FLUTTER SIMULATION INPUTS 







The input vector is; 

f>. 


{•}• 




(C . 38) 


The first term, <5 c , is input to X^q (equation (C.37) 
after being multiplied by the gain (846) ^ in the [b] matrix. 


The noise input, / into the gust terms X 43 through X 


and X^g is multiplied by the constant a 




w 



46 


m 


TT L. 


[B]- 


w 


The output equation is: 

(C. 39) 

The matrices [Hj and [G] contain the modal response 
information: 


{hj = 0] {x} + 0){x} 



m 




h = 

i 

i=l 

rn 

(x,y) 

q i (t) 

(C. 40) 

• 

h = 

t 

i=l 

♦i (x,y) 

q i (t) 

(C. 41) 

h = 

S 

i=l 

4>i (x,y) 

q ± (t) 

(C.42) 


If the wing coordinates for the accelerometer location 
are used then equation (C.42) gives the simulated accelero- 
meter signal. The vertical gust and surface movement were 
also output for analysis purposes. 

C.2 Solution Methods - The transient responses for the 
DAST ARW-1 were generated by a modified version of CONTROL. 
In this computer program the transition matrix and its in- 
tegral are used. These are: 
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=: e [A]t = J +[ A ]t + 2 : [A] 2 t 2 + 3! [A] 3 t 3 


r T n n 
+ n! [A] t 


(C. 43) 


/: 


fA](t T) dx= It + 21 [A] 
[A] ( n " 1 > t n 


t 2 + - [a] 2 t 3 + li 

3 : 


(C. 44) 


The solution of the state equations becomes : 


|x(t)| = 


«Mt) <x f 


lo 


4>(t-x) u*(x)dx (C. 45) 

* * , 
where u (x) = [b] -[u(x)} 

The eigenvalues of the system were in good agreement 

with BWC's results. 

The closed loop responses were generated using a 
simplified FSS control law supplied by BWC. This simplified 
control law matches the complete FSS control law very well 
up to about 400 rad/sec. A block diagram of this symmetric 
control law is shown in figure C.2. 

The accelerations in the individual modes for the 
cases examined in this report are presented in the following: 
Figure C.3 Subcritical Open Loop Responses 
Figure C.4 Subcritical Closed Loop Responses 
Figure C.5 Supercritical Closed Loop Responses 
Figure C.6 Supercritical Closed Loop with Turbulence 
For purposes of comparison to the simulated accelerometer 
signal and control surface movements shown in figure 6.46, the 
following responses are also included: 

Figure C.7 Subcritical Open Loop with Turbulence 
Figure C.8 Subcritical Closed Loop with Turbulence 


320 




SERVO 

ACTUATOR 


FIGURE C . 2 DAST I FSS SIMPLIFIED CONTROL LAW 
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FIGURE C . 3 


SUBCRITICAL OPEN LOOP RESPONSES (PAGE 1 OF 2) 
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FIGURE C . 4 SUBCRITICAL CLOSED LOOP RESPONSES (PAGE 1 OF 2) 
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FIGURE C. 4 
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FIGURE C . 5 SUPERCRITICAL CLOSED LOOP RESPONSES (PAGE 1 OF 2) 
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FIGURE C . 6 SUPERCRITICAL CLOSED LOOP WITH TURBULENCE (PAGE 1 OF 2) 
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FIGURE C . 7 SUBCRITICAL OPEN LOOP WITH TURBULENCE (PAGE 1 OF 3 ) 
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FIGURE C.8 SUBCRITICAL CLOSED LOOP WITH TURBULENCE (PAGE 1 OF 3) 
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